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1.1 


Unit8 Finite-Difference Methods II: Stability 


Set Books 
G. D. Smith, Numerical Solution of Partial Differential Equations (Oxford. 1971). 
H. F. Weinberger. A First Course in Partial Differential Equations (Blaisdell. 1965). 


It is essential to have these books: the course is based on them and will not make 
sense without them. They are referred to in the text as Sand W respectively. 


Unit 8 is based on S: Chapter 3, pages 58 to 72. 


Conventions 


Before working through this text make sure you have read A Guide to the Course: 
. Partial Differential Equations of Applied Mathematics. References to Open University 
courses in mathematics take the form: 


Unit M100 13, Integration IT for the Mathematics Foundation Course. 
Unit M201 23, The Wave Equation for the Linear Mathematics Course. 
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8.0 INTRODUCTION 


In Unit 5, Initial Value Problems we showed how to construct both explicit and implicit 
finite-difference schemes for initial value problems, how to determine their local 
accuracy and how to solve the resulting equations in the implicit case. We demon- 
strated that good local accuracy, measured in terms of the local truncation error, is not 
sufficient to guarantee a successful method, and concluded that for parabolic equations 
the mesh ratio is often an important parameter in this respect. In this unit we look 
more closely at the theory of finite-difference methods applied to initial value problems. 


First we introduce the idea of convergence of a finite-difference scheme. We say that 
a scheme is convergent if, as the mesh spacings are reduced, the finite-difference 
solution gets steadily closer to the true solution of the differential equation either at 
a fixed point, or for all points along the furthest time level under consideration. 
Convergence is our prime objective when using a finite-difference scheme, and in 
practice we might deliberately compute successive approximate solutions, using 
smaller and smaller intervals, until we get results to a required precision. We shall 
see that for this to happen the local truncation error, defined in Unit 5, must tend to 
zero as the mesh spacings approach zero. 


Additionally we need the finite-difference scheme to be stable; that is, small local 
errors of any kind (truncation or rounding errors) should not grow unboundedly 
as the step-by-step computation progresses. Thus convergence depends on both the 
finite-difference scheme and the partial differential equation whereas stability depends 
only on the finite-difference scheme 


In Section 8.2 we shall consider sc me methods for establishing convergence in fairly 
simple cases. These methods are also useful for finding the rate of convergence of a 
scheme. A similar process was followed in the case of ordinary differential equations 
in Unit M201 21, Numerical Solution of Differential Equations. We shall then introduce 
an important theorem, in Section 8.3, which relates convergence to stability and the 
behaviour of the local truncation error, and enables convergence to be established 
more easily and in more general cases. Finally, we concentrate on stability in Section 
8.4. 


Notations often differ in the literature; we shall use the following which, for the most 
part, agree with S. 


Ui; is the true solution U (ih, jk) of the differential equation at the mesh point 
(x, t) = (ih, jk). 
uj is the true (or exact) solution of the finite-difference scheme at (ih, jk). 


ü is the approximate solution of the finite-difference scheme at (ih, jk) obtained 


from a computation. 


hj 


Thus, U; ; contains no error, u; j includes truncation errors and ài; ; has both truncation 
and rounding errors. 
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81 THE BASIC IDEAS OF CONVERGENCE AND STABILITY 


To illustrate the ideas of convergence and stability we shall consider the simple 
explicit scheme given by 


Шыжа — Mig н 21: ig (1) 
k 12 
as applied to the initial-boundary value problem 
e eu 
aU nau xe[0.1], te R*. 
Ct сх" 


U(x,0) = f(x) хед, [0.1]. 
00.0) = ШО, 0) = 0 тед". 


Suppose that we know the solution и at all mesh points up to and including those 
at time level j. We can evaluate an approximation и у, to О; j+1 from Equation (1) by 


шша m tij rea 7 ijt ri) P= ha... N-I, (2) 
Uo je = Uy ei = 0, 

where r = k/h? and Nh = 1. For brevity, we write this as 
п = Luj, (3) 


where uj = (U; js Haj- -+s уст) and L is a (linear) operator RYT! — КУ! which 
can be represented by the matrix 


1-2r r 


r 1—2r 
In Unit 5 (Section 5.5) we defined the local truncation error T; ; for the explicit scheme 
(1) as 
EM Ui iat a Ui; m Visi = 2U,; + U; 


ar Ro 
We also introduced the quantity Т? ; given by 
Thies = КТ, 
mU. — Ung МО 2U;; + Uii (4) 
ie. 
Он = LU; + Thy. (5) 


We now illustrate the assertion in Unit 5 that the local truncation errors can accumul- 


ate as the step-by-step process progresses. We can write Equation (5), for the first 
time step, as 


U, = LU, + Tf. (6) 
Repeating the procedure for the second time step, we get 
U, = LU, + TS 
= L?U, + Т + Tt by Equation (6). 


As usual L? = LeL means “apply the operator L twice”. Repeating the process 
again gives 


U, = LU; + T$ 
= D'U, + DTE + LTS + Tf, 
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and we may deduce the general case, 
U, = L"U, + Y LT (19 = 1), (7) 


del 


by induction. This equation is exact; it relates the true solution of the problem to the 
initial data Ug and the local truncation errors. 


When we solve Equation (1) exactly (without rounding error), we obtain 
u, = Lug, 


u, = Lu, = Lug, 


=[" 
u, = L"ug, 


and since the initial data are given by the true solution at time level zero, that is 
Ug = Ug, we have 


u, = L'U,. (8) 


Thus we can obtain an expression for the global truncation error U, — u, in the 
solution to the finite-difference scheme. From Equations (7) and (8) we obtain, for 
the explicit scheme (1), 


n 


„= Y LU 9) 


jai 


Whether or not these local errors accumulate numerically will depend on the effect 
of the operator L and the nature of the Т;. The local truncation error T, ; may be 
expressed (Unit 5, Section 5.5) in terms of derivatives of the true solution and we can 
infer the nature of the T?;. 


We now concentrate on the rounding errors in the computed solution of Equation (1). 
We assume first that the initial data цо = U; cannot be stored exactly, and so 
the computer begins with an approximation ii; 9 to u; o satisfying 


To = Up + Ro, 


where К, 9 is the initial rounding error. Since there will be a rounding error introduced 
at each stage of the computation, the computer will evaluate uj, , as 


Dj}, = Lū; Ra, (10) 
Thus at the first step we compute 
1, = Lū, + R, = Lu; + LRo + R,, 


and in general 


= Lu, + У OR, 
о 


= LU, + Y L'R. (11) 
=0 


The difference between the computed solution ii and the true solution U of the 
differential equation can be obtained from Equations (7) and (11) and is given by 


= X KTT] = R), (12) 


where for convenience we have set Tj = 0. Thus, in general, both the Tz; and the Rij 
can cause a numerical accumulation 1 in the total error. However, in circumstances 
where the error due to the Tt, does not grow (see, for example, the solution to SAQ 15 
in Unit 5) it is still possible for the error due to the R; j to grow numerically. 


Equation (12) is the expression for the global. or total, error of the explicit scheme (1). 
In practice we cannot calculate the right-hand side of this equation since the 7; ; 


PDE 8.1 


and hence the T£; depend on derivatives of the true solution which we do not know, 
and the R,,, depend almost randomly on the computing device we use. Our best 
bet therefore is to investigate the nature of the operator L to find those conditions 
under which the right-hand side of Equation (12) remains manageable in the sense 
that it is "bounded". Now the essence of stability is that there should be a limit to the 
amount by which any error can be amplified in the numerical procedure. We see that 
the right-hand side of Equation (12) is a sum of error terms Te -Rj operated on by 
some L^. We would hope that the magnification of any error is limited, i.e. that there 
is a real number K such that 


ЕД < КЕД (13) 
for p = 0,1,2,..., and all local errors Ej, where [ЕД = max IE; jl We now say that 


the scheme (1) is stable if there exist Ло and kg such that (13) holds for all p, all local 
errors E; and all h < lig, k < ko. 


Before taking this idea any further we shall introduce the idea of convergence. An 
ideal situation for our numerical scheme would be for the computed solution ii to 
approach the true solution U along some fixed time level as we steadily reduce the 
mesh spacings to zero. Applied to our explicit scheme (1) this would clearly require 
that the sum on the right-hand side of Equation (12) should approach zero as both 
hand k approach zero. Expanding Equation (12) as 
п n 
U,-ü,- Y L'Tt - Y L'R; 
j=0 j=0 

clearly reveals the fact that we can never be assured that ||U, — 9,1 approaches zero 
as h, k ~~ 0 because the second sum depends not only on h and k but also on the com- 
puting device. It is only when К, = 0 for all i,j that convergence in the above sense 
can be obtained. We therefore have to amend our ideas, and we now define a finite- 
difference scheme as convergent when the true solution u of the scheme (without 
rounding errors) tends to the true solution U of the differential equation along some 
fixed time level as we steadily reduce the mesh spacings h and k to zero. Since we 
look at a fixed time level, t say, the fact that л and k tend to zero means that п (= t/k), 
the number of time steps required to reach the fixed time level, increases unboundedly. 
It follows from Equation (9) that to prove convergence we must impose conditions 
on Land TT; so that 


lim у L'T* = 0. (14) 


hk 1 
Now, recall that by definition 
КТ, ТФ, 


ij fje 
for our explicit scheme, so ме want 
anal 
limk У L'77'T;— 0. 
hka jy + 
Now by the generalized triangle inequality we can write 


a жеў 


k Lu EE a-j-YT | 
p j à IL Till (15) 
and if our scheme is stable, we obtain, by Equation (13), the result that 
n-1 amt 
ку LIIT, <k j 
rs J È, КОЁ (16) 


If we let M denote the maximum absolute value of the T; ; then Equation (16) becomes 


n-i 
ky L'UTI| S kn KM. 


* jo 
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However, kn — t and so this inequality becomes 
B1 

k €i 

EY peo, 


i=0 


<КМ. 


Now, if we have the additional condition that 
jim, Ti;20 for all i, j, 


then certainly M also approaches zero. That is, the quantity КМ approaches zero, 
and so the left-hand side of Equation (15) approaches zero and the scheme is con- 
vergent. This latter condition, that the local truncation error T; ; (orkT#,, ,)approaches 
zero as the mesh spacings Л and К tend to zero, is known as consistency (or 
compatibility). We can now say that if our explicit scheme (1) is consistent then 
stability guarantees convergence. This relationship holds for other schemes and 
forms the basis of an important theorem which we shall meet in Section 8.3. 


Although S does not consider local truncation errors in his descriptive treatment 
of stability, we have illustrated both here and in Unit 5 that they play an important 
part in this subject. As a result there is a difference of emphasis between S and what 
has been said so far in the correspondence text. Therefore, we are not asking you to 
read S: page 55, line 1 to page 58, line 13. 
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8.2 CONVERGENCE IN SIMPLE CASES 


8.2.1 The Explicit Formula 


Notes 

(i) S: page 58, line — 13 
e is called the global truncation error in our terminology. 

(i) S: page 59, lines 4 to 12 
These equations are exact using the "remainder" form of Taylor's Theorem which 
states: 


If u is n + 1 times differentiable on the interval [x, x + h] then there is 
a 0 c (0, 1) such that 


n 


h 
u(x + h) = u(x) + hu'(x) +- + a u(x) 


prt 


T ny u”t x 4 Oh). 


(л < 0, then [x, x + h] is replaced by [x + h, x].) 


This form of Taylor’s Theorem is a stronger result than the one we quoted in 
Unit 5. 


In obtaining Equation (3.4) from the equations on lines 5, 6 and 7 in S: page 59, 
use has been made of the Intermediate Value Theorem which states: 


If a function f is differentiable on the interval [а, b], then for any real 
number. y between f'(a) and f'(b) there is at least one point x, e (a, b) 
such that f'(x) = y. 


In particular, there exists x, € (a, b) such that 
P(X) = 3[/ (а) + F'O. 
Hence 
1fo?u eu aU 
Wo Au 0% + 0,h t) + ax ae өл} = ы + ды) 
where 0 < 0, < 1 and 0 < Ө, < 1, for some 0, such that —1 < 0, < 1. 


For further details of these theorems and related topics, sec Appendix 1 of 
Kreider, Kuller, Ostberg and Perkins, An Introduction to Linear Analysis 
(Addison-Wesley, 1966). 


fiii) S: page 59, line — 1 
The finite-difference scheme uses the same (known) initial data as the differential 
equation. Therefore, и; = U; and so е; = 0 for each i. 

General Comment 


We can relate the discussion in S to what we know about local truncation errors. 
We notice that Equation (3.3) on page 59 contains the expression 


Uij+1 Ui + rQU;, — Uii; — Ui y) 


which is T#;,, by definition. In this case Ты = KT, 5. where Т; ; is the local trun- 


cation error of the explicit scheme. Therefore, we may write 
ёра S еъ +L — 20е, + ге, + hE. 


(Compare this with Equation (3.4) in S.) 


10 
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Let M be the maximum value of the modulus of the Т. for all i and j; then 
[ыза] € reci + 0 — len + е] + kM 

which is the expression given in S: page 59, line —8. 

In Unit 5, we showed that 
Tij = O(h?) + О(К) 


and hence T; ; approaches 0 as л, k ~ 0. Therefore, M approaches 0 as л, k ^0, and we 
conclude the convergence proof in a manner similar to that given at the top of 
S: page 60. 


Similar techniques can be used for a simple explicit method applied in a hyperbolic 
case. This is the basis of the following SAQ. 


SAQ 1 

Let U be the solution of the differential equation 
U PU 
AG CAR xeR, г> 0, 


subject to the initial conditions 


U(x,0)= f(x)  xeR 
05,0 = g(x)  xeR 


The differential equation is approximated by the explicit finite-difference scheme 
Wi jer = Mita Маја + ois 


where the mesh ratio p? = k?/h? = 1, and the derivative initial condition is approxi- 
mated by the forward-difference formula 


QUio wu. — Yio, 
дї К 
(а) Show that e; , = O(k?) 
where e; , = U;; — tij- 
(b) Show that 6 541 = езу — еу. + е, 
when U is infinitely differentiable. 


(c) Express e; j+, as a sum of errors on the zeroth and first time levels, and deduce 
that и converges to U as k tends to zero. 


(Solution on p. 23.) 
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8.2.2 The Courant-Friedrichs-Lewy Condition for Hyperbolic 
Equations 


In this section we explore by means of SAQs how the characteristics of a hyperbolic 
partial differential equation can help to determine whether a given finite-difference 
scheme for the equation is satisfactory. Throughout this section we shall refer to the 
following pure initial value problem and finite-difference scheme. 


Problem A 
Solve 
ou eu 
ip as? xeR, г> 0, 
U(x,0) = f(x) xeR, 
ди 


7060) = g(x) хек, 
д 
where c is a positive constant. 
Scheme B 
Up jar — 2ujj + Mij- = AAU уз = Qu; + uii.) ieZ, jeZ*, 
where the mesh ratio p? is given by p = k/h. 
We have seen in Unit 2, Classification and Characteristics that the characteristics 
of Problem A which pass through the point (X, ?), > 0, are 


1 
t—ī= +-(х – х). 
d 


The domain of dependence of the point (X, 1), 1 0, for Problem A is the shaded area 
in the diagram, specified by 


lena e 61-15-59 {жт ш E 
È 


LA 


> 
* 


We define the numerical domain of dependence of a point (x, 7), for a given finite- 
difference scheme, as the set of mesh points whose numerical domains of influence 
(Unit 5, Section 5.2) include (X, ?). Thus the numerical domain of dependence of a 
general mesh point (X, f) = (ih, jk) for Scheme B is given by 


(ih jj 20, j«j-Gü-) jep-ü-1). 


and is shown by crosses in the following diagram. 


12 


SAQ 2 

On what part of the given initial data do the following solutions at (X, Т) depend? 
(a) the solution to Problem A 

(b) the solution to the numerical Scheme B applied to Problem A 


(Solution on p.26.) 


SAQ 3 


Sketch a diagram showing the domain of dependence of a point (X, t) for Problem A. 
On the same diagram sketch the numerical domain of dependence of the same 
point (X, Т) for Scheme B assuming that p > 1/c. What happens to the relative positions 
of the two domains if the mesh ratio p is altered so that p < 1/с? 


Under what conditions are these domains of dependence the same? 


(Solution on p. 26.) 


5404 


If the numerical domain of dependence of a given point (X. Т) lies within the domain 
of dependence of the same point, what happens to the analytical and numerical 
solutions at (X, 7) if changes are made to the initial data at points along the initial 
line which are: 


(i) in both domains of dependence; 
(ii) in the domain of dependence but not in the numerical domain of dependence? 


(Solution on p. 27.) 


540 5 
What general condition must be placed оп the mesh ratio for а finite-difference 


scheme applied to a hyperbolic equation to ensure that all the initial data are used? 
Apply this condition to Scheme B when it is employed to solve Problem A. 


(Solution on p. 27.) 


We have seen that the domain of dependence of a point in the solution domain of the 
difference equation must not lie inside the domain of dependence of the same point 
in the solution domain of the differential equation. We have not, however, investigated 
the case when the numerical domain of dependence of a point contains properly 
the domain of dependence of that point. We now quote a result for this case. It has 
been shown, by Courant, Friedrichs and Lewy. that the contribution to the solution 
of the difference equation at the point (X. ?) from the extra initial data outside the 
domain of dependence (along AB and CD in the next diagram) tends to zero as л, К -«0 
keeping the mesh ratio p? constant, where (X, 7) remains a fixed point. 
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characteristics 
of 
difference 
scheme 


o 
differential 
|| equation 


zu i Eun initial line 
A B С р 


The conclusion is the Courant-Friedrichs-Lewy condition (or just the C.F.L. 
condition) which states that a necessary and sufficient condition for the convergence 
of a difference scheme for a hyperbolic equation is that the numerical domain of 
dependence of a point must include the domain of dependence of the same point. 
SAQ6 
Find the C.F.L. condition for the difference scheme 

шка 205 + щур = p^ (uuu — 20; + ша) 
when applied to the hyperbolic equation 

QU soU 

A cs 

at? ox? 
where c is a positive constant. 


(Solution on p. 27.) 
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8.3 CONSISTENCY AND LAX'S THEOREM 


The local truncation error gives us a measure of the local accuracy of finite-difference 
replacement of a given partial differential equation. We might expect that if we refine 
the mesh, that is, make the mesh spacings h and & smaller, we would reduce the local 
truncation error and would have a more accurate scheme. In particular we would 
like to have 

lim Т; = 0. 

hk-0 
If this condition holds then the finite-difference scheme is consistent (or compatible) 
with the given partial differential equation. Unfortunately, this is not always the case 
unconditionally, as the following SAQ illustrates. 


SAQ7 
Show that Du Fort and Frankel's three-level explicit scheme 


и 


ыж —ULjci — Meg — Шјжа T Шр + Uaj 


2k h* 
has a local truncation error given by 
RU, PU, kU; 
= = +. 
6 а? 12 ax* h? oo? 


when applied to the partial differential equation 
oU — eu 
et Ax? 


and that therefore the scheme is consistent with this differential equation only if 
hand k approach zero in such a way that 


Using the expression for the local truncation error, can you say for what partial 
differential equation the scheme is consistent when the mesh spacings are always 
chosen so that k = ch for a given constant c? 


(Solution on p. 28.) 


5408 
Is the Crank—Nicolson scheme compatible with the partial differential equation 
ди PU, 
a 0^ 


The local truncation error for the Crank-Nicolson scheme is given in SAQ 14 of 
Unit 5. 


(Solution on p. 29.) 


SAQ9 
Is the simple explicit scheme 
u 


ijel = Wiig T Uie quj 7 Uae 


compatible with the partial differential equation 
eU _ eu, 
ar? ax?” 

HINT: See SAQ 15 of Unit 5. 


(Solution on p. 29.) 
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We mentioned in Section 8.1 that there was a relation between stability and conver- 
gence. For the simple explicit scheme applied to the initial value problem 


ыз xeR. 1» 0. 
cet 
U(x. 0) = f(x) xeR. 


we illustrated that compatibility and stability are sufficient to guarantee convergence. 
This condition in fact holds for other problems and forms the basis of the following 
result. 


LAX'S THEOREM 


Suppose we are given a properly posed linear initial value problem and a linear 
finite-difference approximation that is compatible with it. Then the finite-difference 
scheme is convergent if it is stable. 


This theorem applies to both partial and ordinary differential equations. It has been 
proved generally for initial value problems in ordinary differential equations and 
applies in both homogeneous and nonhomogeneous cases. For partial differential 
equations the theorem has been proved only in special cases, although even in cases 
where it has not been proved numerical results seem to indicate that it still applies. 


In the remainder of the unit, therefore, we shall concentrate on stability. This is not 
to say. however. that convergence analysis is of no practical use, since from it we can 
obtain an estimate of the rate of convergence of any particular scheme. 


Note that some authors give a definition of stability less restrictive than ours, for 
which the converse of Lax's Theorem holds, i.e. for compatible schemes convergence 
implies stability. 
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8.4 THE ANALYTICAL TREATMENT OF STABILITY 


8.4.0 Introduction 


Stability analysis deals with the growth of errors in a computation. In the finite- 
difference method we have a recurrence relation approximating a partial differential 
equation and we are interested in whether or not the recurrence relation induces a 
growth of errors as the computation proceeds. Stability analysis of recurrence rela- 
tions is considered in Unit M201 7, Introduction to Numerical Mathematics: Recur- 
rence Relations. For our purposes, however, the recurrence relations involve two 
independent arguments i and j and the analysis is somewhat different. In Unit M201 7 
we were concerned only with the growth of local rounding errors: here we have local 
truncation errors as well. We must determine how and when any local errors 
accumulate. 


In its simplest form stability analysis consists of studying the effect of deliberately 
introducing errors into each stage of a computation and observing what happens 
as the solution to the finite-difference scheme steps forward in time. If the errors 
increase without bound we say that the scheme is unstable. On the other hand if the 
errors remain bounded or even die away then we call the scheme stable. It is often 
the case that a scheme is stable under certain conditions and unstable under others: 
we then say that the scheme is conditionally stable. (Note the analogy with partial 
instability for the corresponding problem in ordinary differential equations discussed 
in Unit M201 21, Numerical Solution of Differential Equations.) 


8.4.1 The Matrix Method 


Notes 


(i) S: page 61, line —5 to page 62, line 3 
We have already stated that there are two forms of error which can accumulate: 
local truncation errors and rounding errors. Truncation errors occur when we 
discretize, that is, replace a derivative by a difference formula (whether in the 
equation or in the initial or boundary conditions). Rounding errors can occur 
at any stage in a computation including the inaccurate storage of initial and 
boundary data. We must, thercfore, include the effects of boundary and initial 
conditions in our stability analyses. 


Our method is to examine the subsequent behaviour of errors at points along 
some time level (not necessarily the initial line) as the scheme progresses in 
time. For convenience we begin our calculations atj = 0. 


(ii) S: page 62, lines —9 to —3 E 
This sentence may be omitted as it merely restates the property of linear inde- 
pendence in another form. 


(iii) S: page 63, lines 4 to 8 
Clearly, 


Aly, = Av, 


v) 


(vi) 
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so that 


= Ale, = AJ oy 
e; = Alea = А У CY. 
s=1 


which is Equation (3.6). 


The result 
N-1 
e= У сй, 
sz1 


shows that the only quantities involving j are powers of the eigenvalues À,. 
If we now introduce the spectral radius of the matrix M, 


p(M) = max |41, 
i 


where 4; are the eigenvalues of M, then the previous result implies 


Nok N-1 


lel € X. led CC = LoCOP У [ШАР 


s=1 


If p(A) < 1 then [p(A)} approaches 0 as j increases. Therefore е; approaches 0, 
the errors die out, and the scheme is stable. If p(A) = 1 then 
N-1 
lel < У, lel Mv, = a constant. 
5=1 
The errors are bounded and the scheme is again stable. However, if p(A) > | 
then |le;]| increases unboundedly with increasing j and the scheme is unstable. 


Thus the criterion for a stable scheme is 

p(A) <1. 
In fact, this criterion is valid even if an eigenvector basis does not exist. 
S: page 64, lines 1 and 2 


This is shown in Unit M201 10, Jordan Normal Form (Section 10.2.2), when f 


is a polynomial function. It is easily shown that if 2 is an eigenvalue of B then 


27} is an eigenvalue of B7}. 


S: page 64, lines 10 and 11 
We have 


r< [а EST з=1,2,...,М—1. 


Now, since 0 < sin?(sz/2N) < 1, the minimum value of 


Jed] 


is 3. Therefore the scheme is stable whenever r < 4. 
S: page 65, line —6 


This equation is derived in the appendix (Section 8.7); the result quoted in S: page 
63 for Ту; is a special case. The notation i(j)k means 


litjhic-2j2gk 


General Comment 


The equation discussed in the reading passage is homogeneous and has homogeneous 
boundary conditions. As a result its finite-difference replacement yields the homo- 
geneous matrix equation 


uj, = Ац. 


For problems which are nonhomogeneous in the differential equation or boundary 
conditions or both, the matrix equations which arise do not have the same simple 
form as above. This is the subject of the following SAQ. 


SAQ 10 
(a) Show that the matrix equations which arise in the solution of the problem 
aU өш | 
= = = О<х<11> 0, 
et Ox 
О(х,0) = 0 0<х&1, 


U(0,r) = f(t) t> 

U(L, £) = g(t) 12 
using the explicit scheme 

Ui jer = Uag +L 2) + rusa 
can be written as 

пу. = Au; + dj 


where the d, are dependent on the given (known) boundary values introduced 
at the jth step. 


(b) If an implicit scheme is used in place of the explicit scheme in (a) the matrix 
equations can be written as 


Buj,, = Au; + dj. 


Show that it is sufficient to investigate the eigenvalues of the matrix B^! A for 
stability if we assume that the known boundary values are given exactly. 


(с) Ifwe use a computer it is likely that the boundary values cannot be stored exactly. 
Determine a criterion for stability if we take into account local errors in the 
boundary values in d;. 


(Solution on p. 29.) 


SAQ 11 
S: page 90, Exercise 2(a). 
(Solution on p. 30.) 
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8.4.2 Two Useful Theorems 


It is quite easy to develop finite-difference schemes and write them in matrix form: 
what is not so easy is to obtain the eigenvalues of the matrices for the stability analysis. 
So far we have met matrices whose eigenvalues have been known explicitly, but 
unfortunately this is not often the case in practice. However, we do not need to know 
all the eigenvalues of a matrix—we require only the eigenvalue with largest modulus 
(the spectral radius of the matrix). The next reading passage introduces two useful 
theorems which can be used in estimating spectral radii. Although it is necessary to 
know only the statements of the theorems a better understanding will be gained by 
studying their proofs (which are not too difficult). 


It should be recalled that if these techniques fail to give suitable estimates for the 
spectral radius of a matrix we can always turn to Unit M201 30, Numerical Solution 
of Eigenvalue Problems to find numerical methods for obtaining eigenvalues; these 
methods are often both fast and accurate. 


Notes 


(i) S: page 66, line —7 
Itis more common to refer to Brauer's Theorem as Gerschgorin’s Circle Theorem. 


(Ш) S: page 68, lines 5 and 6 
Note that h, and h, are physical constants and not mesh spacings. 


(iii) S: page 69, line — 1 
Some logic has been omitted in deriving this expression: the case a,, = 1 — 2r, 
Р, = 2r leading to the condition that r < 4 by Brauers Theorem has been left 
out. The overall stability condition should be 


r< min -— Е ee 
5 2+ hjóx' 2 + hx ^ 2 


However, as л, > 0 and h, > 0, this condition is equivalent to the one on the 
last line of S: page 69. 


SAQ 12 

S: page 91, Exercise 4. 
(Solution on p. 31.) 
SAQ 13 

S: page 92, Exercise 5. 
(Solution in S: page 92.) 
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8.4.3 Von Neumann’s Method 


We have seen (S: page 62. line 7 and SAQ 10) that the formula for propagation of 
any form of error introduced at every mesh point along j = 0 is given by (the homo- 
geneous equation corresponding to) the finite-difference scheme. For example, con- 
sider the explicit scheme 


Us jar = Usa. 7 Qr — Dui; t rusa 

and start the computation with the vector of values ug instead of ug along j = 0 where 
Е, —u, — uj. 

Then 
Epari = ГЕ, — Qr — Ер + TE pig: (1) 


[Note that we have changed the notation for the error vector from e to E in order 
to avoid confusion with the exponential function which we shall be using. We have 
also changed the notation for a mesh point from (i,j) to (p. q) so that we can use i 
later as the conventional symbol for JC 1] 


Stability analysis boils down to solving recurrence relations of the type (1) subject 
to given initial values Е, о (p =0,1,..., N). Von Neumann's method, by analogy 
with the solution of linear partial differential equations with constant coefficients, 
attempts to obtain a solution to the recurrence relation by separation of variables. 
That is, we look for solutions in th: form 


Epa = ХТ, 0) 


where X, represents a function of p only and 7, a function of q only. 


SAQ 14 
By putting E,, = X,T, verify that Equation (1) can be solved by the method of 


pa 
separation of variables. You need not solve the resulting equations for X, and T. 


(Solution on p. 31.) 


We can obtain information about the solution by the method of separation of variables 
when the recurrence relation is linear and has constant coefficients. The method is 
outlined below. 


Consider the initial values E, as values of some error function F, with domain 
[0, 1], at x = ph (p = 0,1, 2,..., №). Suppose F can be expressed as a finite Fourier 
series. As we saw in Unit 6, Fourier Series it is often manipulatively more convenient 
to use a complex Fourier series. We therefore assume that 


N 


Бх) = Y A" хе[0, 1). 
п=0 


Since we аге given the E, о we have the N + 1 equations 
N 
F(ph) = У Ag" = Е о pad; Lesa, (3) 
п=0 

from which the 4, may Бе determined. We are looking for solutions to Equation (1) 
which are of the type (2) and reduce to Equation (3) when q — 0. As the recurrence 
relation is linear we need only look at the behaviour of a single harmonic. Moreover, 
since the A, are constants we shall only be concerned with solutions of the form 


m 
Epa i T. 


where В = пл. If solutions of this form exist then inserting the expression for E,, 
into the finite-difference scheme will yield an ordinary linear recurrence relation with 
constant coefficients for Т,. which will be satisfied by T, = & (Unit M201 7), for 
some choice of č. Hence our quest is for solutions of the form 

E,, = ё?”*ё. (4) 


рч 
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The remaining question is what happens to the term E,, as q increases? Equation 
(4) indicates that Ep Will not increase as q increases provided that |4 < 1, i.e. 

Swe 
This is von Neumann’s criterion for stability. 
We apply this criterion to Equation (1) to get 

oibrheatt — peiPhp+ Dga Qr — 1)erhga 4 yel Lhe 
Division by e#?"&4 gives 

Ё = пей! — (2р — 1) + re Bh 
Recalling that 

ей! + eh = 2 cos Bh, 
we obtain 
= 1 + r(2 cos fh — 2) 


1 — 4r sin? (2. 


Since 0 < sin? (81/2) < 1, von Neumann's criterion is satisfied provided 0 < r < i. 


fan 
| 


Il 


(Compare this result with that obtained by the matrix method, $: page 64.) 


SAQ 15 
S: page 93, Exercise 6(b) and (с). 
(Solution in S: page 94.) 


SAQ 16 (Optional) 
S: page 94, Exercise 7. 
(Solution on p. 32) 


SAQ 17 
S: page 91, Exercise 3(b). 
(Solution in S: page 91.) 


SAQ 18 
Use von Neumann's method to find the stability criterion for the difference scheme 
Upgt1 20,4 + Upg—1 =T (up 44, — 2и, + иа) 
applied to the hyperbolic equation 
QU _ „ди 
en 7 © б” 
where c is a positive constant and r? = k?/h? is the mesh ratio. 
Compare the result with the C.F.L. condition for this problem (SAQ 6). 
(Solution on p. 33.) 
The following important points should be noted concerning von Neumann’s method: 
(a) The method applies only to linear difference equations with constant coefficients. 


(b) Von.Neumann's criterion is always a necessary condition for stability. However, 
it is sufficient only in certain cases; in particular, this is always true of two-level 
schemes. 
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(c) The method applies, in theory, only to initial value problems with periodic 
initial data. 


Notwithstanding the limitations stated above, von Neumann's method, being easy 
to apply. is used in many cases where in theory it need not apply. Often its predictions 
on the stability of a scheme are found to be correct, or at least very nearly so, when 
the scheme is used in practice. 


PDE 8.5 


8.5 SUMMARY 


This unit has concentrated on the theory associated with finite-difference methods 
applied to initial value and initial-boundary value problems. We have discussed and 
illustrated the following points. 


There are two types of error which must be considered in the analysis of finite-differ- 
ence schemes. Local truncation errors arise from the fact that a finite-difference 
scheme is only an approximation to the partial differential equation. We get local 
rounding errors whenever inexact arithmetic is used to solve the finite-difference 
equations. 


It is possible for the local errors to accumulate to such an extent that the global error 
(the difference between the computed and the true solutions) becomes so large that 
the computed solution bears no relation to the true solution. In general we cannot 
measure the local error. and hence the global error, because the local truncation 
errors are expressed in terms of derivatives of the true solution which we do not know, 
and the local rounding errors depend on the computing device used and are essentially 
random quantities. 


For confidence in a computed solution there are two requirements of any finite- 
difference scheme. The first is that a reduction in the size of the mesh spacings should 
produce a more accurate approximation to the true solution of the partial differential 
equation. It follows that in the limit, as the mesh spacings are reduced to zero, the 
computed solution (obtained with exact arithmetic) should tend to the true solution 
of the differential equation. In order that this convergence can take place, a finite- 
difference scheme must also be consistent with the dilferential equation in the sense 
that the local truncation errors tend to zero as the mesh spacings tend to zero. 


The second requirement is that any errors introduced into the numerical process 
should not accumulate or be amplified as the scheme progresses. We illustrated that 
the global error will depend on the nature of the local truncation errors, the local 
rounding errors and also the finite-difference scheme itself. Stability analysis is the 
study of whether local errors will or will not accumulate drastically as the step-by-step 
scheme progresses. We showed two commonly used methods for investigating the 
stability of linear finite-difference schemes; the matrix method and von Neumann's 
method. Lax's Theorem tells us that, for properly posed linear initial value problems, 
a linear finite-difference scheme which is both stable and consistent is also convergent. 
Even though Lax's Theorem tells us that we need to investigate only the stability 
and consistency of a scheme, convergence analyses are not redundant since they 
are the only means of estimating convergence rates. We therefore studied ways of 
establishing convergence in the case of two simple explicit methods. We also used 
the fact that hyperbolic equations have two real characteristics, the properties of 
which gave rise to the Courant- Friedrichs- Lewy condition for convergence of 
finite-difference schemes for hyperbolic equations. 
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8.6 SOLUTIONS TO SELF-ASSESSMENT QUESTIONS 


Solution to SAQ 1 


(a) 


(b) 


We have 


glih) = CU p, 0) а 0 UMP, 
Ct 
and so 
ша = fit Kg. 


where f; represents f (ih) апа g; represents g(ik). Taylor's Theorem gives 


eU. 5 
Uia = Ui + k= + 000) 
Е 
= fi + kg; + О(?). 
Thus. 
ёа = Uia ша = О(К?). 
In fact, 


le; 4l S 3M. 
where M is the maximum value of |&* U/0t?| for x € К. t e (0. А]. 
Substituting и; = Шу — е, into the finite-dilTerence equation 
Шыжа = Hisig T Mage FMR 
we obtain 
biga = Creag ®ы- tO ug Giga — Gin yt Ша Uu ü) 
By definition the local truncation error T; ; is given by 


Ui; -2U;; + U Visi Чы + Uii 


Ts ij ыта | 
EJ k? h? 
or 
kT j= Uij+i „> Шай + Uij-ı Ui-ij. 
since = k, and so Equation (1) becomes 
Ciji = Craig T бы-1 t Cic + ЮТ, (2) 


In SAQ 15 of Unit 5 we showed that for this problem 7; ; = 0 for p = k/h = 1 
when U is infinitely differentiable. Hence, from Equation (2), we get 


ejs = Crag T бы-1 + е1. (3) 
Replacing i by i + 1 andj + 1 by j in Equation (3), we get 

Cnty = Cisaj-i — Fist joa + быз. 
and similarly, replacing i by i — 1 and j + 1 by j in Equation (3). we get 

Crag 7 €5j-1 — €i-1j-2 + Ci-ij-iv- 
Substitution of the last two equations into Equation (3) yields 

Cigar = Creager T Cij-i + Ciji-3j-1 T Cieij-2 T Ci-1j-2 (4) 


and we have expressed e; j+; їп terms of errors along the j — 1 and j — 2 levels. 
We can write down expressions for the errors along thej — 1 level in Equation (4), 
by employing the recurrence relation (3), to express e;;,, in terms of errors 
along the j — 2 and j — 3 levels as 


Cigar = Cieaj-a + Cie gj-2 t Cicij-2 + Ci-3j-2 7 Cie 24-3 


з. (5) 


— €jj-3 — @-» 


ы 
D 
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We can see that Equation (3) gives us e; j+, in terms of two errors along the jth level 
and one error along the (j — 1)th, Equation (4) gives it in terms of three errors along 
the (j — 1)th level and two along the (j — 2)th, and Equation (5) gives it as four 
errors along the (j — 2)th level and three along the (j — 3)th. In general we can see 
that we can write e; j}; in terms of n + 1 errors along the (j — n + 1)th level and n 
errors along the (j — n)th. We could prove this by induction. When n — j we obtain 
the relation 


€, +1 = (/ + 1 errors along level 1) — (jerrors along level 0). 


Since the same initial values along j = 0 are used for both the differential equation 
and the difference equation we have е; о = 0 for all i and, from the result of part (a). 


we have 
le; i] < $^ M 


where M is the largest value of |22 U/0t?| in the first time interval. Hence, 


less € G+ DIEM 
or 
le; € 4jk?M. 


Now let jk = t, some fixed time; then |e, | < $tkM, which approaches 0 as k ~ 0, 
provided 0^U/Ot? is bounded in the first time interval. 


Therefore the scheme is convergent. 


Solution to SAQ 2 


(a) The solution to Problem A at (X, 7) depends on the data on that part of the initial 
line included in the domain of dependence of (х, i). 


(b) The solution to Scheme B at (x, ?) depends on the data at those mesh points at 
level 0 which are included in the numerical domain of dependence of (x, i). 


Solution to SAQ 3 


characteristics numerical 
of Problem A characteristics 
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The domains of dependence (for Problem A) are shaded in the diagrams; the hatched 
areas are the numerical domains of dependence (for Scheme B). 


When p > 1/с the numerical domain of dependence of the point (x, f) lies within the 
domain of dependence of the same point. When р < 1/с the numerical domain of 
dependence lies outside the domain of dependence. The numerical domain of depend- 
ence varies as the mesh ratio p is altered, and in general the disposition of the numerical 
domain of dependence of a point relative to the domain of dependence of the same 
point depends on whether p < l/c or p > I/c. 


The domains of dependence of (х, 7) for Problem A and Scheme B coincide when 
p= le 
Solution to SAQ 4 


In the diagram the initial line segment BC corresponds to (i) and the initial line 
segments AB and CD correspond to (ii). 


characteristic 


A B C D x 


Any changes to the initial data along BC will cause changes in the analytical solution 
and in the numerical solution at the point (X, 7). However, changes along AB or CD 
will cause a change in the analytical solution but not in the numerical solution. 
Therefore, no matter how good an approximation the numerical scheme is to the 
hyperbolic equation, the solution of the numerical scheme cannot respond to changes 
in the initial data in the way that the analytical solution does. 


We have shown the characteristics of the hyperbolic equation as straight lines since, 
in Problem A, c is a constant. The argument of the solution to this SAQ still holds 
when the characteristics of the given hyperbolic equation are curved. 


Solution to SAQ 5 


By the result of SAQ 4 the mesh ratio for a finite-difference scheme applied to a 
hyperbolic equation should be chosen so that the numerical domain of dependence 
includes the domain of dependence. 


For Problem A, since the characteristics are straight lines, the mesh ratio should be 
chosen so that 


1 
ae 
BS 
Solution to SAQ 6 
The characteristics of 
QU PU 
at 


ax? 


through the point (X, 7) are given by 


tp alg on. (1) 
с 
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The numerical characteristics of the scheme 
thar — DM p.m р? (шару Hij + шед 
through (X. 7) are given by 
ї—-ї= p(x- X) (2) 


The C.F.L. condition requires that the numerical domain of dependence of (X.T) 
should include the domain of dependence of the same point. This implies that the 
slopes of (2) be less than or equal to the slopes of (1) in magnitude. That is 


p 


I^ 
ate 


Solution to SAQ 7 
Using Taylor's Theorem. we find that 


ей, КОО, 


Шуу — Uij-1 = Mo т 53 + 0089), 
QU,, k* e*U, М 
Us Uii 92U, IO + IL =i 4+ O(k9), 
ijti + Uij-a ij t дг? + 12 art + 0n) 
апі 
А Е , ht 8*U,, 
Ui, i, + Uiig = 2U,; +h EM "b 12 LPS + O(n’). 


Therefore, the local truncation error is given by 


Ui jet n Uii = Uiij;— Ui jay a U,i-i + U; 
2k h? 


=j 


Tij 


"m k PUn 
ij _ ij + 6/52 
D ae то oae tOO + 00m) 
ёш eus. kUJ 12 eu, k? eu, 
et 6x? 6 2 12 ё т д? 


+ O(k*) + O(h*) + O(Kk*/I3). 
Then, for the differential equation 
ди cU 


ё ёх? 


the local truncation error is 


к? BU; h? 9*U;; к? Ж 


Te" FF тобе tg ag + OK) OU) + Ouchih?). 


We can see that the expression for the local truncation error cont 


Hila ains a term involving 
/h^, which does not approach zero as h, k ~ 0 unless 


0. 


lim = = 
y h 
If k/h = c then the equation for the local truncation error becomes 


ij 


U.. OU. ppm 2 да 2 
20 = PU " k? OU, | e Uj, 2 9^U,, 
д бх* 6 ar 12 ах ? 


+ O(Kk*) + O$) + Olek?) 
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which approaches 0 as Л, k ~ 0 if U satisfies the differential equation 
eU eu eu 


Qr e ёг? 


which is hyperbolic! Thus the scheme is consistent with this partial differential equa- 
tion provided the mesh spacing is always chosen such that k = ch. where c is constant. 


Solution to SAQ 8 


From SAQ 14 of Unit 5 the local truncation error of the Crank- Nicolson implicit 
method applied to 


-eU PU 
à ox 


is given by 


Т; = O(k?) + O(h’). 


d 
from which we see immediately that T; approaches 0 as k, h ~ 0 and the scheme is 
compatible with the differential equation. 
Solution to SAQ 9 
From SAQ 15 of Unit 5 the local truncation error of the simple explicit scheme 
Ui jer = Mica, Up quj — Маја 
applied to the differential equation 
eu eu 
is given by 
Tj = 0, 


and the scheme is compatible with the differential equation. 


Solution to SAQ 10 


(a) Thejth matrix equation is 


Ujer (1— 2r) r un 1, 
Ua jad r (1-2) r из; 0 
= + 
| Uy- isi r (1—2r) | usi rg; 


which can be written as 
и; = Au; + dj. 

(b) The Crank-Nicolson formula is given by the equation on the last line of S: page 
64. The inclusion of boundary values means that we have to add the column 
vector (rf;, 0. 0..... rgj to the right-hand side of this equation which now 
takes the form 

Buj,, = Au; + dj. 
(The same form would be obtained if we used some other implicit method.) 
This last equation can be written as 
-- pL -1 
иы = B^!Auj + B dj 


provided det B # 0 (which is the same as saying that it is possible to solve the 
given matrix equation). 


(c) 
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If we now introduce errors at every point along j — 0 and begin the computation 
with the vector uj instead of ug we see that for = 0. 1. 2... 


fac ВА + Bod, (1) 


provided that there are no errors in the d;. 
Putting 
canc 

ej = uj uj * 
we obtain 

еы = В Ae; 
which implies that 
= (B^! Aeg 
which does not involve any boundary values, and we need only consider the 
eigenvalues of B^! A. 


i 


If we do introduce errors into the boundary values we have to write Equation (1) 
as 

uf,- B Auf + Bo 'd* 
where the d? are the boundary values d; with errors introduced. Putting 
e;-uj— uf once more, we obtain 

еа = B^! Ae; + B^ (dj — df). 
Putting b; for B^'!(d; — df), we see that 

4-1 


e; = (B! Aeg + Y. (B^ Ay !7*b,. 
k=0 


Again, the error vector is bounded and the scheme is stable provided the spectral 
radius of B~'A is less than or equal to one. 


Solution to SAQ 11 


The difference equation may be written as 


(1 + 2ra)ju; рар rats jti H Mie je) = s 


ij 
which, in matrix form, is 
1 + 2ra —ra LET uy dij 
—ra 1+2ra —ra [M LN 0 
—ra | 4 2ға —ra Ux 2,41 Uy-2j 0 
—ra 1 + 2ra usc jei Wye dy, 


where the d;; are obtained from the known boundary values. We can write this 
equation as 


U — raTy- uj, = uj + d; 


where Ty_, is a square matrix of order N — 1 given by 


30 


-2 1 
1-2 1 
227 
Tac [3d 
1-2 1 
i 29 
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Thus 
ча = (I — raTy a) tu, 0 — гаТу-) 14, 


provided det(/ ~raTy_,) #0. The eigenvalues of (J — raTy .,)^! are (using the 
eigenvalues of Ty_, given in S: page 63) 


1 


Ty s=1,2,..., N—-1, 
NEN! 
1 + 4rasin ES 


and these are strictly less than 1 in modulus for all positive values of r, so that the 
difference scheme is unconditionally stable. 


Solution to SAQ 12 

The scheme of part (a) yields 
пз = Au, 

where 


1-2r r 


r 1-2r 


and the solution follows by the method of S: page 92. 


Solution to SAQ 13 


See the solution in S: page 92. 


Solution to SAQ 14 
Putting E,, = X,T, in 
Еа+1 = TE nei — Qr — DE,, + PE азах 
we obtain 
Ха = Х.Т, 0r — ХТ, + rX, AT, 
Dividing by X 
T, 


per Xa —Or— DX, РХ 


T, X 


4 Р 


ply. We obtain 


and. since the left-hand side of the equation is a function of q only and the right-hand 
side is a function of p only. the equation is separable as 


Tai 
rN.) — r — DX, Хр = EN, 

where č is a constant. (In principle we can solve these two recurrence relations by the 

methods of Unit M201 7.1 

Solution to SAQ 15 


See the solution in S: page 94. 
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Solution to SAQ 16 


Since the error function E satisfies the original difference equation (with homo- 
geneous boundary conditions), substitution of E,, = ейт gives 


ейюга+1 ойт = n (gift Diga _ өйте 4. uil heal 4 kbelbhe, 
п 
Division by e?" gives 
ka, _. Я 
ё-1 = m (e^ — 2 + e) + kb. 
2 


Hence, 


4k [ 
E=kb+1- тзн" (07. 


For stability we require |£| < 1 and therefore we need 


4ka . , {Bh 
-1&kbtl— quom B $1. 
Recalling that В = пл and Nh = 1 we require the following. D 
4k 
(i +1 05211 «1  D1Xn&N-i 


= B = 
п ° O2N 
since the boundary values are known. This is equivalent to 


4a пт 
X — sin! —. 
b S ртіп yy 


which we may replace by 


4a л 
h? € —sin? —, 
b 2N 
since 
nna 
sin? — < sin? — <п<М—1. 
ШУ; sin” oy І £n€&N-—1 
Hence, 


as required. 


4ka . , nx 


i) -1<k l- -= — 
(ii) b 91 aw 


is equivalent to 


„пи А 
4ka іл? sN S (kb + 2)h^, 


i.e., : 
2a sin? < k + ih. 
NN ko? 
so that 
P sud PERO? aga 
у  2asin ay Hob. or 
k 1 
IO PIS UT DEF l<n<N-1. 
2a sin? — ib 
2N 
Since 
nt 
1 > sin? zy En&N-l, 
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we have 


as required. 


Solution to SAQ 17 
See the solution in S: page 91. 


Solution to SAQ 18 


Putting Ep = ет" in 
Eye) ~ 2Epa + Epa- = FC, tig = Epa + Epaia) 
we obtain 


eifiPhgas 1 на 2e'Prhga + ei nhe 1 m r?c? (gpt Опра _ 2gilirhga + ейі Dhgay 


Непсе, 
ё 2+ EH сей —2 4 еті), 
ie. 
. a [f 
6-2 + ёё! = —r?cl4sin? (5. 
Thus 


8-23 |: — 27202 sin? (=) č+1=0, 
which gives 


či č = 1 — 27202 sin? 


Now, for stability we require 
ё &1 1= 1,2. 


If г2с2 sin?(Bh/2) < 1, we have 


"n 132 {Bh oV [fh 2o a [PRAE 
= elg + ERE] 
Či. Ča = 1 — 2r^c? sin |5 + 2ire sin 5 1 — с? sin 5 


апа |é,| = 1 124,2. 
On the other hand. if rc? sin? (fhy2) > l we get [| > 1. 
Therefore 

lél <1 1= 1,2, 


holds only when r?c? sin?(Bh/2) < 1, and so is true for all f if r?c? < 1. Thus, von Neu- 
mann's condition for stability is k/h < 1/c. From SAQ 6 the C.F.L. condition for 
convergence is also k/h € 1/c. 
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87 APPENDIX 
The Eigenvalues of a Common Tridiagonal Matrix 


The following is a derivation of the formula given in S: page 65 for the eigenvalues 
of the tridiagonal matrix of order N given by 


a b 
c a b 
c a b 


c a 
where b and c are real and have the same sign. 


By the definition of an eigenvalue we have the matrix equation 


a b ДД vy 
c a b ГА ПА 
= 4 + 
c a Ъ [оуу Dy 1 
c а Uy Vy 


where (r,..... vy) is the eigenvector associated with the eigenvalue 7. We see that 
this matrix equation is equivalent to the set of equations 


av, + bv = Avy, 
cv, + av, + bo... = Av, r=2,3,...,N—1, 
соус + aly = Avy. 
If we define vg = гу; = 0, the above equations can be written concisely as 
CU, + av, + bou = Av, r= ТО, N, (1) 
which is a second-order recurrence relation subject to the boundary conditions 
Uo = tye = 0. (2) 
The solution of (1) is given by (see Unit M201 7) 
v, = dip, + dps 
where p, and p; are the roots of the quadratic equation 
bp? + (а – jp - c-0 


and d, and d; are arbitrary constants. Now, using a trick from Unit M20/ 25 (Section 


25.4. we make the substitution д — a = 2(bc)? cos. So the quadratic equation 
becomes 


+ 
: c c 
ma X et = 
p 15 cos 0 p- т = 0, 


and the roots are 


Py. Pz = Hi (cos Ü + isin 0). 
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The general solution to the recurrence relation (1) is therefore, by De Moivre's 
Theorem, ш 


"T [р © (dy cos 0 + d, sin r0) 


where d, = d, + d, and d, = i(d, — dj). 
The boundary conditions (2) yield 


43 20 
and 
c (N+ 1)2 
B d4 sin(N + 1)0 = 0, 
for which 
5л 
—MÁ s=1,2,...,N. 
d N+1 ? 


Thus the sth eigenvector has its rth component given by 


c. ors 
QV = b sin В 


d*. зл fe Е 257 d E Мл 
n =s] a |=] si р 
ПЕЧИ ON ST b N+! 


with corresponding eigenvalue 2, given by 


that is, 


Àj-a + 2(be}* cos т e 
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It is essential to have these books: the course is based on them and will not make 
sense without them. They are referred to in the text as S and W respectively. 


Unit 9 is based on W: Chapter V, Sections 27 and 28. 


Conventions 


Before working through this text make sure you have read A Guide to the Course: 
Partial Differential Equations of Applied Mathematics. References to Open University 
courses in mathematics take the form: 


Unit M 100 13, Integration 11 for the Mathematics Foundation Course, 


Unit M201 23, The Wave Equation for the Linear Mathematics Course. 
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This book contains an excellent account of Green's functions for boundary value 
problems in Chapter 12 (Sections 12-9 to 12-11). The approach is from a slightly 
different viewpoint from that adopted in this unit, in a style which will be familiar 
to students who have taken M201. 


9.0 INTRODUCTION 


The present unit treats ordinary differential equations, subject to initial conditions 
and boundary conditions, from a fresh point of view. Firstly we summarize some of 
the principal properties already known to you. 


The general linear equation of the second order is of the form 


аи du 
a(x) dd + a(x) а= + acx)u = f(x), (09; 
and the associated homogeneous equation is 
2y dv 
alx) qs + аул + ах) = 0. (2) 


The solution space of Equation (2) has dimension 2; that is to say, any three solutions 
of (2) are linearly dependent, and there exist two solutions which are linearly inde- 
pendent. If we call a given linearly independent pair of solutions v, and v;, then for 
any constants c, and сз. 


CU, + 6202 
is a solution of (2). Moreover, every solution of (2) is of this form (Unit M201 9, Differ- 
ential Equations 11). 
Let н, be any solution of (1). Then (1) has the solution set 

{ити = up + сүрү + C202; C1, C2 ER]. (3) 


Practical methods are needed for obtaining up, v, and t3. When аз, a,, ag are con- 
stants, v, and v, are easy to obtain (Unit M201 9), but otherwise. a special investigation 
is normally needed. To find a particular solution u, we may employ intelligent 
guesswork for simple cases, or use a systematic method such as Variation of Para- 
meters (Unit M201 11, Differential Equations IT). 


One solution in the set (3) may be singled out if initial conditions are supplied, of the 
form 


ulg) = ug, u'(x) = uy. 


This requirement says that the solution curve sought must pass through the point 
X = a, и = ug, and have slope equal to u, at x = о. The Existence and Uniqueness 
Theorem (Unit M201 33, Existence and Uniqueness Theorem for Differential Equations) 
states that, under quite general conditions on аз, а, ag and f (which are satisfied 
for all the specific equations considered in this unit), there is one, and only one, 
solution to Equation (1) which satisfies the initial conditions. 


In this unit we shall consider principally another class of problems, using the above 
theory as a starting point. These are boundary value problems—to find a function 
u satisfying 


аи” + au + dou = f in (а, P), 
with 
u(x) = а, u(f)- b. 


We have met such problems in Unit M201 25, Boundary Value Problems. Here we 
require, if possible, a solution curve passing through the points (2. a), (В. b). 


The Existence Theorem does not say anything directly about these problems, there 
may be no solution, there may be one and only one. or there may be an infinite 
number. Simple examples of each possibility are the following: 


(д ui tu=0 in (0, я), 
u(0) = 0, u(x) = 1. 


This has no solution. 


(i) uw +u=1 in (0, л). 
и(0) = 0. u(x) = 0. 


This has the unique solution u(x) = 1 — sin x — cos x. 


(ii) «tus in (0. x). 
и(0) = 1. u(x)-l. 
This has an infinite number of solutions, u(x) = 1 + C sin х, for any value of C. 


The occurrence of these possibilities does not really depend on delicate properties 
of f and the coefficients; the boundary value problem has completely new features 
differing from those of the initial value problem. The method of Green's functions 
which is discussed in this unit enables these possibilities to be investigated, and allows 
the solution to be written down when it exists. The method applies also to differential 
equations in domains of dimension greater than one (that is to say, partial differential 
equations) and Unit 10, Green's Functions I will deal with this topic. 


Our method of attack is as follows. We begin by constructing a formula for the solution 
to the nonhomogeneous initial value problem, in terms of an influence function. We 
then use this result to give us some leverage on the nonhomogencous boundary value 
problem, for which we obtain a formula in terms of a Green's function. The Green's 
function is essentially a continuous, but not smooth, "solution"' to the associated 
homogeneous problem. The basic Green's function solution is refined (in Section 
9.2.2) to deal with the most general type of (unmixed) boundary conditions, and a 
physical interpretation is given in Section 9.2.3. 
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9.1 INITIAL VALUE PROBLEMS 
" 


Notes 


(i) W: page 117, first paragraph 
The transformation leading to Equation (27.1) is obtained as follows. When the 
term (pu^) is expanded, the equation becomes 


pu" + ри + ди = Г іп (0, В), 

and we require this equation to Бе equivalent to the original equation 
au” + Би + cu = Е in (0, ff), 

that is to say, to have the same solution set. 


Assume that a and p have nonzero values throughout (о, f/); then we require the 


equations 
Ё 4 
и" +? ME af 
p р р 
апі 

PISTE F 
Wt cd pa Se 
a a a 


to be equivalent on (о, f). This is achieved if the coefficients are the same, i.e. if 
Тоша FoF € 2X4 1 
p a p a р а (0 


The first of these is a differential equation, with a solution 


=E) 
p(x) = exp fl RE 


ls 


is any indefinite integral of b/a. (The one chosen in W is 


where 


ra dE " 


x alc) 


In the footnote it is explained that when lim ač) is zero it might be necessary to 
ma 

choose another lower limit—in other words another indefinite integral.) 

The functions д and f are now obtained from Equations (1). 


(п) W: page 118, line —5 
We have К = — р(х) М [v, , v2](x) for all x e (ж, f) where W denotes the Wronskian 
(Unit M201 9). Clearly p(x) is nonzero throughout (а, $), and the Wronskian is 
nonzero since v, and v, are linearly independent. Hence, К + 0. 


(iii) W: page 119, line 6 
In particular. К is independent of any particular choice of v, and v3. 


(iv) W: page 119, lines —11, —10 
The function 
wx | RESSO 


is the particular solution of the equation 


(риу + ди = f 


PDE 9.1 
which satisfies the (homogeneous) initial conditions u,(%) = 1102) = 0. The 
solution set of the differential equation is given by 

{ии = uy + сүгү + C303 1C, C2 E RJ. 


where гү and гу are any two linearly independent solutions of the associated 
homogeneous equation 


(pr) + qv = 0. 


The system 


(риу +qu=f — in (a, P) 
ula) = а, ula) = Б 


has а unique solution for every pair of constants а, b (Unit M201 33), and this 
solution can be expressed in the form и, + сус, + сз: for some pair of constants 
€. C2. 


We find the constants from the following equations : 
ula) = (a) + cv (x) + catala), 


ie. 

a = cv, (x) + csvs(o), (*) 
and 

u'(a) = ula) + eril) + cavale), 
Le. 


> 
Ш 


cywila) + cavala). (*) 


General Comment 


The derivation of the influence function given here is similar to the method of 
Variation of Parameters described in Unit M201 11. In the notation of that unit, the 
equation on W: page 118, line 4 is essentially of the form 


VAX) = yiexeiQ) + уз(х)сз (х) 


where foreknowledge has provided the author with the correct form of the functions 
c, and сз. The following table illustrates the relationship of the notation used in W 
and that in Unit M201 11. 


M201 Ww Й 
y+ ayy + ау = (риу + 40 = f 
ау, Gg. h р'/р, 410, f /p 
Yi Y2 01, 02 
У, w/K 
X C А 
Green's function K(x, г) influence function R(x,&) 
for an initial value 
problem. 
SAQ 1 : 


Reduce the following equation to self-adjoint form: 
(1 + x*y"(x) + 2xu'(x) + Au(x) = 0 -1<х<1 
(Solution on р. 24.) 
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5402 


Reduce the following equation to self-adjoint form: 


и"(х) + 1 u(x)- (x) 2 F(x) 0Ocxc«l. 


х 


(Solution оп р. 24.) 

5403 

Find the influence function for (ће equation 
u"(x) — u(x) = е“ x > 0, 

and obtain the solution satisfying 
и(0) = 1, w(0)20. 


(Solution on р. 24.) 


5404 
Reduce the equation 

x^w'(x) — 2xu(x) + 2u(x) = xt x»1 
to self-adjoint form. 


Determine two linearly independent solutions v, and v, of.the associated homo- 
geneous equation, and confirm that р(х) [0 (х)оз(х) — v2(x)v,(x)] is constant. 


Obtain the influence function R, and use it to solve the system 
x^w'(x) — 2xu'(x) + 2u(x) = х“ x1 
u(1) = 1, u'(1)- 0. 
HINT: To solve the homogeneous equation, look for a solution of the form x — x". 


(Solution on p. 25.) 
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9.2 BOUNDARY VALUE PROBLEMS 
9.2.1 Homogeneous Boundary Conditions 


In Section 9.1 we saw how the initial value problem: 
(риу + qu=f і (a, P), 
u(a) = 0, w(a)= 0, (1) 


can be solved by constructing the influence function R(x, 2). This requires only a 
«nowledge of any two linearly independent solutions, v, and р, of the associated 
homogeneous equation (рг) + qv = 0. The solution of (1) can then be written in the 
orm 


u(x) = | R(x. ë (2) 
f the initial conditions in (1) are nonhomogeneous, so that и(о) = а, u'(x) = b where 
a and b arc not both zero, the problem can be solved by adding suitable multiples 
of v, and v; to (2), as іп W: page 119. 


n this course we are concerned as much with boundary value problems, where a 
condition is imposed at each end of an interval, as with initial value problems. A 
typical example is 


(риу +qu=f і (0, р), 
u(x) = 0, (6) = 0. (3) 
This system has homogeneous boundary conditions, but the equation is nonhomo- 


geneous. We shall find again that the solution can be represented as an integral, this 
time in the form 


В 


u(x) = | (х, č) f Ode, 


a 


where the function G turns out to be constructed from solutions of the homogeneous 
equation, as the influence function was. 


This method provides a useful way of displaying the solution, which, as we shall see 
later on, has a strong physical motivation. It also enables the question of the existence 
and uniqueness of solutions of the system (3), and of the analogous system with non- 
homogeneous boundary conditions, to be resolved. Remember that the existence 
and uniqueness of solutions of ordinary differential equations has so far been con- 
sidered for initial value problems only. 


Notes 


(i) W: page 121, line 11 
The case D — 0 is discussed later. 

(0) W: page 121, line —8 
You will be asked to fill in this manipulation in SAQ 5. 

(їп) W: page 122, lines 4 to 6 
It is very important to grasp the meaning of the statement “G is symmetric". 
Merely to exchange x and č in the algebraic expressions in Equation (28.2) 
(W: page 121) does not return G(x, 2) unless they are also exchanged in the 
inequalities. The following notation is helpful. Let 


1 
G*(x,č) = Kp Ovi) = vay G9] [0 i(x)e3(B) — va(x)0,(B)] 


and 


1 
G (xé)- Kp Level) — vaxo Q0] [o Eel) — vy )e, 09] 


(iv) 


PDF 9.2.1 
for (x. £) e [a fj] x [a, 3]. Then 


es Joxe € 
бв ihi н 


AN IK 
Fa 


If we exchange < and x in the whole expression above (so that 2 < x becomes 
x < С and vice versa), we get 


E ё, 
сх. 
Ви 
б, х) = G^(x.£) and G^(&.x) = Сх, E. 
ѕо that 
3 G^ (x.&) xS, 
БЕ м E& x. 


which is the same as G(x, ©). 


W: page 122, Equations (28.4) 
Note that here "d/dx" really means “0/0х, č being kept constant"; č has the 


status of a fixed parameter in these equations, which are satisfied by the- function 
of one variable 


хк» G(x, ё) x € [о, fj]. 


These equations will from now on be adopted as an alternative definition of G. 
They serve also as a practical way of obtaining G without remembering (if you 
ever could) the formula (28.2). The conditions (a) to (d) should be remembered; 
they follow directly from the definition of G, as you will be asked to show in 
SAQ 6. It can be shown that these equations are sufficient to determine G. if it 
exists. In particular. this implies that G is independent of any particular choice 
of v, and г. 

W: page 122, Example 

We give a picture of the Green's function for this case. It is clearly its own mirror- 


image in the plane x = ё 


G(x, à) 
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SAQ 5 


Complete the “algebraic manipulation” required to establish W: page 121, lines 
—7. —6. 


(Solution on p. 26.) 


SAQ 6 


Confirm that for each č e (x, В). С. defined by Equation (28.2). satisfies the conditions 
(28.4) in W: page 122. 


(Solution on p. 26.) 


SAQ7 
Consider the system 
u+u=0 іл (0, 3л). 
(0) = ия) = 0. 


Show that there is no solution except the trivial solution и = 0. Find the Green's 
function for the system and confirm its symmetry. 


(Solution on p. 27.) 
SAQS 
Consider the system 
и'+и=0 in (0, x), 
u(0) = u(x) = 0. 


Show that there is a solution, and that no Green's function is defined by the boundary 
value problem (28.4) in W: page 122. Check also that D = 0 for the system. 


(Solution on p. 28.) 
5409 
(i) Let 
D(f.g) = f(3g(9) — f (B)g(o) 


for any two functions fand gon (a, B]. Let (v, , v2} and {v3, v4) be any two pairs 
of linearly independent solutions of the equation 


(pv) + qu = 0 in (o, 8). 
Prove that 
D(v,, 03) = 0<> (гу, v4) = 0, 


so thal the conditions D = 0 and D + 0(W: page 121, line 11) are independent 
of the choice of solutions used. 


(ii) Prove that if the system 
(ри) + qu = 0, 
u(a) = u(B) = 0, 
has a nontrivial solution, then D = 0. 
(11) Prove that if the system 
(ри) + qu = 0, 
ща) = u(fl) = 0, 
has no solution (except for u = 0), then D + 0, 
(iv) Deduce the converses of (ii) and (iii), 
(Solution on p. 28.) 
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It follows from SAQ 9 that the construction of the Green’s function G is possible 
for a boundary value problem if and only if the associated homogeneous problem 
has no nontrivial solution. 


9.2.2 Nonhomogeneous Boundary Conditions 


In the next reading passage we see how to tackle nonhomogeneous boundary con- 
ditions. We then have a brief discussion which enables us to define a Green's function 
where the boundar nditions i involve d. ives. Finally, an example is treated. 


es 
(i) W: page 123, lines 5 and 10 
The notations 


Са 


ёс 
áp x) and ag ^ p) 


are ambiguous, and if care is not taken wrong results may follow. We illustrate 
the correct process as follows. 


We have 
0G* 
Be OE) = ев) — о) n] 
"m x [0 ,(x)v2(B) — vGv B] x> é, 
(x С) = Es 
` * (x,¢) = PIN) = 02(х)01(0)) 


x [vi(C)va(B) — vée) x«é 
We now define the function of one variable which maps x to 
eG 


=z (x,a) = ae 9 хе [о, В] 
д Kex 
| 0G*, , 
дЕ (x, € Ж 
= ЕБ [w(a)vala) — volalo] l)o) — va(x)v,(B)] 
= ij 6900 = бг] 
from the definition of K. Therefore 
2 (о, x) = — M 
ge OD pla) 
and 
ёа 
E (B. a) = 0. 
For 
3G 
25 (В. В) and ae x, В) 
дё 


you must be careful to use G^ rather than G*. Concisely, 
~ 


(ii) 


(iii) 


2G АС” А act 
с n= дын (x. 2) Е (B, %) = z— (х. 2) 

CS GG yea ee х=й 
^ 2G- eG 2G- 

ia. py = Ep ĉc ip py = oS (x. fl 

* ‹ EM e сс х=й. 


W: page 123. lines 11 to 16 А 
We want to use the known solution of the nonhomogeneous equation with 
homogeneous boundary conditions (28.1) in РУ: page 120 to solve the completely 
nonhomogeneous system (28.5). The method used depends on the fact that 


eG(x,é 
Uu iX 0—3 oe ) x e[a. ff] 
ФЕ |а 
satisfies 
(py — qu, = 0, 
1 
по) = — #0, u,(f) = 0, 
pla) 
and that 
xe[x, f] 
satisfies 
(рш) — qu, = 0, 
1 
@)=0, (0) = — = + 0. 
аты 


Since u, and и; are linearly independent, all solutions of the differential equation 
are of the form 


u = up + Au, + Bu, 


where 
f 
u(x) -[ GG, 0 f (dé 


is a particular solution. In fact, the particular solution u, satisfies the homo- 
geneous boundary conditions, which makes it easy to adjust the constants 
A and B so that u(a) = a, (В) = b. 


W: page 124, Equations (28.6) 


The special cases obtained by making some of the constants jt, 413,0, Огоз zero 
should be noted. If gu, = ji = 0 we have Dirichlet conditions; if c, = a, = 0, 
Neumann conditions. 


SAQ 10 


Show that the system 


u + 4и = 1 in (0, л), 
и(0) = 1, u(x) = 0 


has no solution. 


(Solution on p. 29.) 
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SAQ II 
Solve the system 
w- u= -f in (0, 1), 
u(0) = 0, w'(1) = 1. 
by obtaining first the appropriate Green's function. 


(Solution on p. 29.) 
SAQ 12 


Confirm that, in the notation of W: page 123, dines 4 to 10, the Green's function of the 
example in W: page 122 satisfies 


eG e 
fogli Sao =o, 
ЕЕ ec 

6G 


== (0, 1) = 0, буш —1. 
eg 


Use this Green’s function to write down the solution of the system 
и = -1 їп (0, 1), 

u(0) = u(1) = 1, 
and check your answer. 
(Solution on p. 30.) 
You have seen in the examples and SAQs that, when constructing a Green’s function 
in practice, it is more convenient to use Equations (28.4) rather than the explicit 
formula (28.2). It is of interest to see how the symmetry of the Green’s function follows 
directly from the properties stated in Equations (28.4). (In practice it is often more 
efficient to use symmetry instead of continuity—(28.4c)—as in the examples in 
W: pages 122 and 125.) We shall also show directly that the Green’s function deter- 
mined by Equations (28.4) yields the solution of the nonhomogeneous problem 
(28.5). The advantage of this treatment is seen when we consider the more general 
system (28.6), for which Weinberger does not give the explicit form of the Green's 


function. Again, in the next unit, when we come to consider Green's functions for 
partial differential equations there is a close analogy with Equations (28.4). 


We first prove a simple lemma which you have seen previously in Unit M20/ 25, 
Boundary Value Problems. 


LAGRANGE'S IDENTITY 
Let L be the operator defined by 
L:u > (риу + qu 


where u is twice continuously differentiable on some interval. Then 


[en = rLu) = рш — си). 


Proof 
We have 
foe = rlu) = || [и(р]' = cpu] by definition of L 
= [ир — ёри) — [ie = ри] integrating by parts 


= puc — ru’). 


You may recognize that we have essentially used this result in showing that p(x) x the 
Wronskian is constant over the interval (W: page 118) and in obtaining the criterion 
in W: page 124, line 6. We now use Lagrange's Identity to obtain two further results. 
The methods of proof are somewhat involved. and we are not suggesting that you 
should be able to reproduce them "blindfold". However, we hope that you will 
appreciate the ways in which we use the defining properties (28.4) for the Green's 
function in the midst of the manipulations. 


We shall first prove symmetry for the function G defined by Equations (28.4): we 
shall show that 


G(x, x3) = С(х;, Ху), 

whenever x, and x, e (a. fj). 

Proof 

Let x, and x; be fixed points іп [æ, 3) with x, < x, and set 
gix > G(x, xj) x € (х, В]. 
By 2X H+ G(x, x3) x € [a, f]. 


We apply Lagrange’s Identity to g, and g, over the interval [«, fj): 


xi х: а 
| (81182 — gabe | {81185 — #2181) + f (8178; — 82181) 


x70 х2 
= С = аз) * [meses = z 


à ху+0 


о B 


+ Е = z 


where we have split the interval to account for the jump discontinuities in the integrand. 
By Equation (28.4a) in W: page 122, Lg, = Lg; = 0 over each interval and the left- 
hand side of this expression vanishes. In addition we note that 


gile) = gx) = 5100) = 8:8) = 0 бот (28.4b), 


gı and g, are continuous at x, and x; 


хэ+0 


gi is continuous at x; from (28.4a, c). 
gh is continuous at x, 
Thus, collecting terms, 
0 = р(х,)а›(х)[81(Х + 0) – gix — 0) 
+ plx2)ei(x2)Lgo(x2 — 0) — eoi; + 0] 
Now, by (28.4d) 


gx + 0) - gl 0) = — 1 
p(x) 
fori = 1,2. Hence, 
82(X1) = gi(xz) 
or 
G(x1, x3) = GG x). 
as required. 
We shall now take the nonhomogeneous boundary value problem, 
(pé) + qu=—f іп (a, f) 
u(a) = а, 
u(B) = b, 


and show that its solution is given by И/: page 123, line — 16 when С is defined by 


Equations (28.4) іп W: page 122, assuming that both (28.4) and the above problem 
have a solution. 
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We set, for some fixed Xe (o, fj]. 
gx G(x, X) — xe[o fi], 
and apply Lagrange's Identity to u and g over the interval (о, fj): 


x Li $-0 
Í (ulg — gLu)+ | (Lg — glu) = [ч = en] eee = el! 


2 FEST 


The left-hand side reduces to 


[ gf. 


since Lg = 0 except at х and Lu = —f over the interval. To evaluate the right-hand 
side we use the conditions 

g(a) = g(f) = 0 from (28.4b), 

g is continuous at ¥ from (28.4c), 


и and и are continuous at X 


и(о) =a, u(f)- b. 


Thus, 
? " 
gf = – р(ођав (о) + p(X)utX)[gX — 0) — g'( + 0)] + p(B)bg'(8). 
By (28.4d) 
gi — 0) - gE + 0) = — 
P(X) 
so that 


p 
щ®) = fer + ap(a)g'ta) ~ bp(B)g'(B). 


а 


Now, G is symmetric, as we have already shown. Therefore 
g(&) = G(&, X) = GG. č), 


and so 


Hence 
р Р 
u(X) = Í G(X ё) f (dé + ара) 52 ® а) — bol) G. B. 
for all X e [«, f]. 


SAQ 13 

By proceeding as on W: page 121, show that a unique solution to the general non- 
homogeneous system (28.6) does not exist unless D # 0, where D is defined two lines 
below Equations (28.6). 


x х a 
nan: (f но, да) = Hox [ КҮН 
dx ls, Ox 


хо 


(Solution on p. 31.) 


SAQ 14 (Optional) 


By proceeding as in the preceding text, show that the Green’s function G defined by 
the system (28.7) is symmetric. 


(Solution on p. 32.) 


SAQ 15 (Optional) 
By proceeding as in the preceding text. prove that (subject to D = О and the existence 
of a solution). the solution to the system 


-j in (2. fj]. 


li 


(pu) + qu 
u(x)- а. WR = b, 
is given by 
ГА 
u(x) = | G(x. (2142 — plaaGtx. з) + pf) b G(x. f). 
where G is the Green's function satisfying (28.7) with ш = —1. e, = 0. wo = 1, 
o,=0. 


(Solution on p. 32.) 


9.2.3 An Intuitive Interpretation of the Green's Function 


The solution of the system 
(риу = ди = — in (e fl). a) 
u(x) = u(f) = 0. 
can be represented (when it exists) in the form 
n 
us) [ ec are. (2) 
А 
By analysing this integral we shall demonstrate that such a representation is to be 
expected as a consequence of the linearity of the differential operator 
L:u > (ри + qu 


and of the linear and homogeneous nature of the boundary conditions. We shall also 
give a direct interpretation of the Green's function G. The argument is not rigorous, 
and there is no need to learn it. 


Partition the interval (2, f] intoa large number of small intervals (2;.. , , £j], 1 1,2,... n, 
where 


y-2565€«à«àxzv x6 - f. 


It will do no harm to assume that the intervals are all equal in length: 


G) 
` 
We consider now the system (1) with the function f replaced by F;, defined as 
" Ше Či <Ë Së, 
Fé) = 4^ : Е 
e n elsewhere. (4) 


A graph of F; is shown in the next figure. 


PDE 9.2.3 


* [T B € 


Let the solution to Equation (1) for this case be denoted by и. Then Equation (2) 
takes the form 


p : 
но) = | GG Ра = GG ESEA О 


Comparing (5) with (3), we see that 


u(x) = У u(x). (6) 


i=. 


Since 
fo) = У Fix), (7) 
ї=1 


Equation (6) states that the solution of (1) is the sum of the solutions of several Sub- 
problems with f replaced by the function F;. f is built up from the F; (i = 1,2,...,л) 
by (7), as shown. the width of the columns being A. 


2 ~ 


> 
a B x 


This result is predictable from the linearity of Land the homogeneity of the boundary 
conditions in (1). In fact, since 


Lu; = (ри) — qu; = —F;. 
u;(x) = uff) = 0, 


fori 2 L.... n, it follows that 


Y Lu; 


\ 


У [pu — ди] 
i=1 


=- У Е = -f from Equation (7) 


with 


n n 


Y ula) = у uf) = 0. 


i=] i= 
But the first equation can be written as 


de MESI] e > 


=1 izl 


so that the solution of (1) satisfies (6). The validity of the argument is subject to the 
supposition that, by taking А small enough. the error in the solution of the system 
when using the step function approximation to f. instead of f itself. can be made 
arbitrarily small. 

We now find the meaning of G(x. ). We see, by applying Equation (5) to the function 
J given by 


(8) 


elsewhere. 


that G(x. &) is the solution to the problem of an input function which is large in 
(&—1. ©] and zero elsewhere, and which has the property that 


В 
[fede =. (9) 


In fact f resembles the delta function of Unit M201 12, Linear Functionals, or the 
impulse function of Unit M201 19, Laplace Transforms. Such a concentrated input 
often has a physical interpretation, as in the following example. 


The system 
u(x) = — EET O<x<l, (10) 


u(0) = u(1) = 0, 


describes the steady displacement n(x) of a string stretched between x = 0 and 
x = 1, under a force per unit length equal to g(x), where T is the tension in the 
string. (Equation (10) can be obtained from the vibrating string problem in W: page 36, 
by assuming there is no time variation.) 


The function g may represent any kind of force distribution. We shall suppose that it 
represents a concentrated force such as might be produced by hanging a weight 
from the string at the point x = č, as shown in the figure. 


vx 


force 


u(x) 


There is no way of producing a force concentrated wholly at a single point—there is 
always a spread of force about the point x = & but physical intuition leads us to 
expect that the exact distribution does not matter much, on the scale on which we are 
viewing the phenomenon, and that the important quantity is the magnitude of the 
resultant force, given by 


1 
| g(x)dx. 
° 


Then the argument above suggests that u(x) is equal to G(x, ë) where G is the 
Green's function for (10) (see Equations (28.4) and the example in W: page 122). 
It is obvious from the figure that the function и shares some of the properties of G: 
except at x = С, we have u"(x) = 0, u is continuous at $ and has a discontinuity in the 
slope there. 


The interpretation of the Green's function as the effect ofa concentrated input function 


to the system permits a physical interpretation to be made of the symmetry of the 
Green’s function. In the case of the string, the symmetry property, 
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PDE9.23 


С(хь, хо) = G(xo. xp). 


states that the displacement at P due to a point force at Q is equal to the displacement 
at Q due to an equal point force at P (see figure). This reciprocity law is of wide import- 
ance in mathematical physics. 


үх 
o 
o 

vx 


0 9 Р I 
| 
1 
1 
1 
1 
1 
1 


Y G(x, хо) 


SAQ 16 


The differential equation 


mx + Кх = 


describes the motion of a particle of mass m attached to a spring of stiffness К and 
acted on by a force f(t) at time г. Find the solution for t > 0 when f(t) reduces to a 
unit impulse delivered at t = 0, the system being quiescent for t « 0. (A unit impulse 
increases the momentum of the particle (= mass x velocity) by one unit. The 
notations X and X are conventionally used to denote the first and second derivatives 
of x with respect to t.) 


(Solution on p. 33.) 


PDE 9.3 


9.3 SUMMARY 
We have seen how ordinary dillerential equations may be expressed in self-adjoint 
form as 
Ци) =f 
where L is the operator defined by 
Щн] = (pu'Y + qu. 


for suitable functions p and q. We have exhibited solutions for this equation with 
initial and boundary conditions. 


The system with homogeneous initial values 
Ци] =f in (0. x). 
ulg) = (в) = 0. 


has the solution 
us = [ etd. 


where R is the influence function or one-sided Green's function. R is given by 


Е vides) — eov (0) 
R(x, č) = ME 2 PARE 
PAX) (ey (x ох) (х) 
where v, and v, are any linearly independent solutions of the associated homogeneous 
equation 


Lie] = 0. 


We saw that R is independent of the choice of v, and v3, and that for any particular 
pair гу, v the denominator is constant throughout the interval. 


The system with nonhomogeneous initial values 
Ци) =f in (0, 2). 
ulg) = ug. 
u(a) = ty, 


has the solution 
чә = [RIESE + сура) eto. 


where the constants e, and сз can be determined by substituting the initial conditions. 


For boundary value problems we must consider the associated homogeneous equation 
with homogeneous boundary conditions: 


Ци) = 0 in (a, ff), 
— gi u'(x) + oula) = 0, (1) 
Haw (B) + aul h) = 0. 


This may have a nontrivial solution v. say. In this case cr is a solution for every 
constant c, and there are no other solutions. There is no Green's function (SAQ 9) 
and the solution of the nonhomogeneous problem does not necessarily exist; if it 
does exist it is not unique. 

If there is no solution to the associated homogeneous problem then there is a Green's 
Junction G, and for each £ є (2, fl) the function x H+ G(x, $) satisfies the following 
equations. 


(а) L[G]-0 х= 


6, 


а 
= fa = 


+ 0,6 
To dx 


х= 


= 0, 


хер 


х= +0 
(c) [5] = 0, 
х=&-0 


- m О EUM 
dx |eaz-0 р(ё) 

С is symmetric: 

G(x. ё) = G(E. x). 
In the case that дү = р = 0. G is given explicitly by Equation (28.2) in W. 
The general boundary value problem 

Ци) = —f — in(o.f), 
= pta) + aula) = a. 
H2u'(B) + ou(B) = b. 


in which neither д, = 6, = О nor jt; = о; = 0, has the solution: 


u(x) = | G(x, ё) Г(&)4ёЁ + са a) + — i bG(x, В) 
fy 
when и, # 0, po #0: 
D А 
= | Gi, fide + aS uy PD, C eg) 
Ж e, бё b 
when 4, = рә = 0; 
Li 6G 
по) = [ awara + PO асрх,о — 00 578 x, py 
a Hy O2 
when u, #0, р; = 0; 
в 0G 
u(x) = | Сбх, ё) f (€)d£ + 28), E (х,а) + me S BGs .p 
x 1 i 


when д, = 0, i; #0. 


кә 
D 


9.4 SOLUTIONS TO SELF-ASSESSMENT QUESTIONS 


Solution to SAQ 1 
This can be done on sight*: the self-adjoint form is 


afi + Teoh aao -l<x<l. 
dx dx 
Solution to SAQ 2 


The choice 


р(х) = e| | | 


will not work, because the integral does not converge (the trouble can be spotted 
from the fact that formally we obtain exp(In x — In 0) and In 0 is undefined). However, 
any indefinite integral of 1/x will do as the exponent; so we choose the simplest. 


[dx 


p(x) = exp | =) = exp(In x) = x. 


We now write, in the notation of W: page 117, 


р(х)с(х) п? 
g(x) = =-— 
a(x X 
рбх)Е(х) | 
fo) a ^ xF(X). 


The self-adjoint form is then 


L аы) — T a(x) = xF(x). 


dx V dx ^ x 


Solution to SAQ 3 


The equation may be written 


ue \ то) = u(x) = е“, 


dx \ dx 
so it is already in self-adjoint form with p = 1.4 = – 1, f(x) = œ. 


Two linearly independent solutions of the corresponding homogeneous equation 
are v(x) = ех, v(x) = е^, We have 


К = p(x)[vi(x)va(x) — vo(x)o(x)] = 2. 


(The fact that K is constant and nonzero is a partial confirmation that the manipula- 
tions so far are right.) 


Then. from Equation (27.4) in W: page 118 
R(x, ё) = 4(e%e7 = — SeT). 


The solution is therefore 
1[* 
u(x) = Hi (етет — ee") dé суе + cse? 
o 


= хе + зет — de") euet е7“. 
The conditions u(0) = 1, u'(0) = O givec, = c, = 4. The solution is therefore given by 


u(x) = 4xe* + de* + Зет“, 


* Assuming, of course, that you have good sight. 
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Solution to SAQ 4 
We want to find p, q, f so that the equation 
x?u"(x) — 2xw'(x) + 2u(x) = xt х> 1 
is equivalent to the equation 
(ри) + qu = pu" + pu + qu = f. 
This is achieved if 


p) O 2 qu) 2 fe) x 
;- 


рх) x p x 


7 


р(х) OX 


The first of these equations is satisfied by 


2 
р(х) = a| f- an| = 
x E 
| = 2 ay 
x 


could be any indefinite integral of — 2/x, and we have naturally chosen the simplest 
one. Then 


where 


g(x) = 2/x* and f(x) = 1. 


The self-adjoint form is therefore 


df dui m 2 1 
db zuel, 
dx \x? dx w 
and it has the same solution set as the original equation № x > 1. 


We require two linearly independent solutions of the associated homogeneous 
equation, and we look for a solution of the form x —» х". Substituting v(x) = x" into 


x?v"(x) — 2xv'(x) + 2e(x) = 0 х> 1, 
we obtain 
n(n = 1)x" — 2nx" + 2x" = 0 xl. 
If there is a value of п satisfying this equation we must have 
0 = n(n — 1) – 2л + 2 
= (п — 2)(n – 1). 


that is to say n = 2 or n = 1. and so, the homogeneous differential equation has the 
linearly independent pair of solutions 


гбх) = x, eax) = х2. 


Then. we obtain 
poor on со бх) — esee] = x7 pg = 23] = - 1, 
which is constant. Equation (27.4) in W: page 118 gives the influence function 


2 _ Ex?), 


R(x. ё) = — (xë 
(We can confirm that for a fixed 2 

Rye = 0. 
and 
$ I 


= Op 


dR 


dX zz 


= ( 


Alternatively, we could have used the homogeneous equation with these conditions 
to determine R.) 


to 
D 
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The solution to the given system has the form 


px 


ux = | (= xë? + dx1H dé + сурх + сух? 
1 


¿=x 
U al HeX + сох? 


£-1 
= (= ixt t dx) -(-dx dM eux cxt. 


After applying the initial conditions a(l) = 1, u(1) = 0. we obtain 


lse 60. 
0=c + 203: 
therefore c, = 2, су = — 1. and finally the solution is 


u(x) = ix = 3x3 + х. 


Solution to SAQ 5 
The expression which we have to manipulate is 


v2(@)0y(x) — ty (a)va(x) 


R(x.& + Б 


К(В, č). 


Since 


гу(х)02(6) — рох) (E) 


К(х,&) = K 


and 
D = v,(a)v,(f) — v;(2)o,(8) = KR, B), 
where K is a constant (W: page 118), this becomes 


DK {220 а)о (0) — v2(a)v,(B)) (0.0)0:(8) — р(х) (@)] 
+ [valv (х) — vi) eax) [2 (B)v2(€) — 0208)0,(2)) 
pk" 1(9)е1(х)о2(8)02(&) + vi (B) vs (E) oso) vo (х) 


= 010)04(6)02(0)02(B) — vi(e)v,(B)u)o4(2]. 


It can be confirmed that this is equal to 
1 К Р 
= xp Quit) = vx(&) v.e] [v(x — va) (8). 
as required. 


Solution to SAQ 6 


We write (28.2) in the form 


G*(x, ё) <x< 
G(x,2) = Ue 9 Ex B. 
G^(x,£) а<х< 6. 
(а) Капа Dare constants, assumed not to be zero. Equations (28.2) may be arranged 


to read 
G(x, ё) = A* (Dux) + В" (8) (х) хе (ё, fl 
G(x, ё) = A" (ё) (х) + B (vj) хе[о, 2], 


where A*, B*, A^, B^ are independent of x. Since v, and v, satisfy (28.42), so do 
G* and G^. 


(b) Direct inspection of (28.2) shows that 
(8,8) = 0, G- (a, č) = 0. 
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(c) Inspection of (28.2) shows that G^ (č, č) = С (č, 2). This condition states that 
G is continuous. 


(d) We have 

dG* _ 1 

dx = gp PELA — 0268) 010901 Los eA — 050) 01(0)). 
Therefore 

dG | dG* 

@х|к= +о С dx k= 

1 
= kp ivi) — vivi (o) Do (29) — v2.01. 

Similarly, 

dG ас” 

dx x-£-0 dx х=. 


1 
= кр! Тр 006) (о)1[о,(8)02(8) — valé) (AN. 
By expanding these expressions апа subtracting we find that 


dG 
dx 


dG 
cero 4х 


х=&-0 

1 
= p 1910999000200) + v. (v5 leal) — vi Oye Qvi va) 
= 01(8)02(8)01(8)0(9)] 


= 51-59% + vy (So; (8)] (0100) 0208) — v (B)v2()] 


But the first bracket equals — К/р(2) and the second equals D, so the result 
follows. 
Solution to SAQ 7 
The general solution is 
и = Asin x + B cos x. 
The condition u(0) = 0 gives В = 0, and u(4m) = A sin }n = A is zero only if A = 0 
too. 


To construct the Green’s function it is better not to copy from (28.2), but to build it 
up from (28.4) in the following way. 


G^(x č) = C(č) sin x 
satisfies G7 (0, 2) = О and the differential equation. 
G*(x, č) = D(é) cos x 


satisfies G* (4n, č) = 0 and the differential equation. 
Continuity at x = € requires 
D(é) cos č = C(£) sin č. 
The discontinuity in the derivative (since p = 1) is given by 
—D(é)siné — C(é)cos— = —1. 
Therefore C(é) = cos č, D(z) = sin ё and 


Glx.2) sin © cos x 0<ёх<х 
х,с) = enm 
(6 cosé sin x 0<х< 


To prove symmetry. we interchange x and é (see note (iii) of Section 9.2.1), and obtain 
К sin x cos ë 0О<х<ё< іл, 

осу = RE 0<ё<х «in. 

which is equal to G(x, č). 

Solution to SAQ 8 

A continuous solution is u = sin x (or any multiple of this). 

The Green's function must satisfy the following conditions: 

dG 

dx? 


(b) G(0, č) = G(x, ё) = 0, 
(c) ер. = 0, 
5-0 


dG 2+0 
= --L 
n Fe 


(The notation in (c) and (d) should be self-explanatory.) 


(a) +G=0 Xe. 


Every solution of (a) on (0, č) which is zero at x = O is of the form A (č) sin x, and every 
solution on (Ё, л) which is zero at x = wis of the form В(ё) sin x. Condition (c) requires 
that A = B, since sinë + 0 У e(0, л), and therefore (d) cannot be satisfied. 


Putting v,(x) = sin x and v;(x) = cos x, we obtain 
D = 01(0)0,(я) — e;(0)o,(x) = 0. 


Solution to SAQ 9 
(i) Any three solutions of the homogeneous equation must be linearly dependent, so 
ру = av, + bus, 
by = cv, + dv; 
for some a, b, c, d. Therefore 
D(v3, рд) = (ad — be)D(v,, v2). 


If ad — be = 0, the equations above imply that v3, v4 are linearly dependent, 
contrary to hypothesis. Therefore ad — bc + 0 and the result follows. 


(ii) Let u, be a solution of the given homogeneous system. Then D(u,, u3) = 0 for 
any choice of из, since u(x) = u,(fl) = 0 by definition. It follows from (i) that 
D — 0 whatever solutions are used. 


(iii) Let u}. u, be two linearly independent solutions of the homogeneous equation, 
satisfying respectively 


u(x) = 0, 18) = 0. 
Then w,(B) + 0, и.(о) # 0, or the given condition would be violated. Now 
D(u, , из) = u(x)us(B) — wi(f)us(o) 
= —u,(B)uz(«) # 0. 
Therefore, by (i), D # 0 for every choice of solutions. 


(iv) The converses state that if, for some choice of solutions of the homogeneous 
equation, D = 0, then there exists a solution of the homogeneous system—for 
if there were no such solution, D 4 0 by (iii); and if D + 0, there is no solution 
of the homogeneous system—for if there were a solution, D = 0 by (й). 
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Solution to SAQ 10 
The associated homogeneous problem 
vU + 4 = 0 in (0, л), 
v(0) = (л) = 0, 


has the nontrivial solution v: x к—> sin 2x. Hence the given system will have solutions 


only if the criterion in W: page 124, line 6 is satisfied. For our problem the criterion 
requires 


T" 
25 || sin 2x dx, 


o 


which is untrue. 


Solution to SAQ 1 


The Green’s function G must satisfy the Equations (28.7) which, for our problem, 
reduce to Р 


pe 76=0 хя=& (1) 
dG 

put B 

Glasto = G|.-;-o = 0 (3) 

юрар i 

dx |x=z+0 dx eo ` m 


A linearly independent pair of solutions of (1) is [e*, e^ *). The differential equation and 
conditions (2) are satisfied by 


G(x, ё) = AQ)e te ree) x>g, 
G(x, č) = BEWE = е7") x SE 


Equation (3) requires that 


А(&) (е1 + е5) — BS — e~’) = 0, (5) 
and Equation (4) gives 
Allee? — e7) BE + е7 = =I. (6 
Equations (5) and (6) give 
б —e 5-1 1-4 
AQ I. В@==——— 


e+e!) ere” 


and the Green’s function is given by 


7 (е — е7) (ех! + е'-*) к 
Mi xc 717. "7 


me 


_ (67! + ete — e) " 


С (х, ё) = Yer ej xx. 


This can now be substituted in the appropriate version of the last paragraph of 
W: page 124. However, it is most unlikely that you will remember such a formula, 
and the following procedure is better in this case. 


Let 
1 
g(x) = | Сүх, ё) 70), 
0 


so that 


g-g--f in (0. i). 
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Since 
dG 


G -— 
kso Xl, 


=0 


for every č є (0. 1) (by definition), g must satisfy 
g(0) = g'(1) = 0. 
Therefore we have to add to g some function h which satisfies the conditions 
h-h-0 1п(0,1), 
М0) = 0, (1) = 1. 


The function x — С (x, 1) satisfies 


dos 
DE (x, 1) — G7(x, 1) = 0, 
G^ (0,1) = 0, 
and 
26- 
a D = 
Therefore, 
ES Like. Ro 
h(x) = Gx) = ete? 
and 
1 nae foe 
-f GG. )/(@4& + т 
- оооло + [ eco area +=. 
ete 


where С? and G^ are given by Equation (7). 


Solution to SAQ 12 
We have p(x) = 1, and 


G* (x, d = č(1 — x) х2 Ё, 
G 
ma (е х9 = 0-х x«& 
Now, 
a 
209-29 o ; 
z х=о 
we therefore consider 
oct 
áp o9 =1l-x, 
and letting x = 0 this becomes 1, which is equal to 1/p(0). 
Similarly, 
6G 
a th 021-120 
Next, we require 
3G ôG- 
will 1 
at (LD = 3 ae (x, 1) ae 


30 


3 


Therefore, we consider 


Os 1) = =x, 


which gives us 


1 
1,1) = -1 = –- ==, 
к ! p(t) 
and 
6G 
ac © 1) = 0. 


Using the formula (W: page 123, line 16), we obtain 


x 1 
ues) = [ a= eae + f x(1— 4 41.1.0 —x)  1.1.(—x) 


= (1 — хр? + x[(1 х) — 301 — х®)] +1 
= -4x5 + 5x41. 
Check: и"(х) = —1,u(0) = l,u(1) = -$4+341=1. 
Solution to SAQ 13 


The function 


- алда 


is one solution of the equation (pu') + qu = —f, though it does not necessarily 
satisfy the given boundary conditions. Every solution is therefore of the form 


_ | 2088 - ьу(х)р(@ 


и(х) 


K f (dE + сүр(х) + Сорох), 


where we have written R(x, £) out in full as on W: page 118, and we wish to choose 
су, €; So that the boundary conditions of (28.6) are satisfied. 


Using the differentiation formula in the hint, we have 


(у= р(х)гу(х) — = x69) ec - [ v(x) 02(6) аы ләде 
` + суох) + суру(х) 
-f vi(x) val) = val) ec yit + iod + сэр). 
Therefore, 


u(a) = c,vi(x) + cavala), 
ula) = civi (a) + carla), 


| vy(f)ux(&) — v(e ë) 


K f (d£ + celh) + сур), 


u(B) = — 


у, ) EE n), j 
wip) = - [ О — (600 C) yee + ci cans. 


a 


The boundary conditions therefore become, after some rearrangement: 
c,[- ur (a) + oe y(o)] + cal- vs) + с11(9)] = а, 
сунот) + озо (В) + 2820500) + 0222000) 
Lba [ (uar (f) + oat (ox) — [urat ) + 03209] n) 


K f Gd. 
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PDE9SAQ 13-15 
These equations have a unique solution с,. сз provided the determinant D of the 
coefficients is not zero. If 
р =0, 


a solution may exist. but is not then unique. 


Solution to SAQ 14 

Let x, and x, be fixed points іп (о. В] with x, <x, and set 
gix — G(x. x4) xefa fj]. 
gj:X — G(x, x3) x € (a. fj]. 


We proceed as in Section 9.2.2, noting that the boundary conditions are now, using 
(28.75), 


— Hy Bila) + agila) = 2808) + ogh) = 0 
fori 2 1. 2. We obtain 
0 = —p(a)[g (algr(a) — ga(o)g1(0] + р(х) (х1) [gi (x + 0) — gi(x — 0)] 
+ р(х2)81(х2)[85(х — 0) — ga(x + 0)] + paf Ces(/g5(8) — gag (] 


and when we "plug in" the boundary conditions the first and last brackets vanish. 
Finally, the jump condition (28.7d) gives 


Оч + 0) – g(x; 0) = – Ee 
for i = 1,2;so that 
82(X1) = 81 (Xp) 
or 
G(x, , x3) = G(x, ху), 
and we have proved symmetry. 
Solution to SAQ 15 
For some fixed х e [a, В] we set 
&:Хҥ—› G(x, X) x e [o. В], 


and apply Lagrange's Identity as in Section 9.2.2. noting that the boundary con- 
ditions are now 


g(x) = g(f) = 0 
u(x) = a, (В) = b. 
Thus we obtain 
f 
f af = р(о)в(ођа + р(Х)щ(%)[е'(Х — 0) — g'& + 0] — р(8)2(8)Ь. 


By (28.7d) 

g& - 0) -gG 40 
р(х) 
so that 


f 
u(x) = [е — р(о)ав(о) + p(B)bg(f). 


Now G is symmetric, as shown in the solution to SAQ 14. Therefore 
g(€) = G(é, x) = G(X, 2). 
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Hence 


b 
(х) = | GR, E) 708) — р(@@)аб(х, о) + p(f)bG(s, В), 
for all X e [a, fj]. 
Solution to SAQ 16 


The velocity x immediately after the unit impulse is delivered is given by 
mx(0 + 0) = 1, 


and the displacement is still zero, х(0 + 0) = 0. Subsequently the displacement is 
governed by the system 


m + Kx =0 in R*, 
1 (1) 
х(0) = 0, 0) = — 
т 


The general solution of the equation is 


к\+ 

“| |, 

т 
and the initial conditions give 


Н 
B=0 and 45) =]: 


KM 
Asin (5 t + Bcos 
m 


m 


Mp 


The system defining R on W: page 119 is analogous to the system (1). This question 
elucidates the remark preceding the first example on W: page 119. 
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PDE 10.0 


10.0 INTRODUCTION 


We begin this unit by looking at Laplace's equation in a circle, solutions for which 
we shall require in the sections following. We go on to extend the ideas of Unit 9, 
Green's Functions I, where Green's functions were introduced lor ordinary dif- 
ferential equations, and solve nonhomogeneous partial differential equations in two 


dimensions. Our aim is achieved, in the first instance, by introducing the finite Fourier 
transform. 


We then try a different approach using ó-functions in order that we may relate our 
results to the definition of the one-dimensional Green's function in Unit 9 and use 
analogous definitions in two and three dimensions. The Green's function treatment is 
used finally to treat problems with nonhomogeneous boundary conditions, using 
Poisson's equation in two dimensions to illustrate the techniques. 


PDE 10.1.1/10.1.2 


10.1 THE FINITE FOURIER TRANSFORM 


10.1.1  Laplace's Equation in a Circle 


Notes 

(i) W: page 100, line —8 
In order to conform with the usage in W, we shall uselogfor the natural logarithm 
function from now on. 


(ii) W: page 100, line —5 
A singular point of the ordinary differential equation 
ау dy 


а(х) тз + pm 


is a point x where either a(x) — 0, or b(x) or c(x) is infinite. Thus, a normal dif- 
ferential equation has no singular points. 


(iii). W: page 101, lines 8 to 10 
Uniform convergence and term-by-term differentiation were discussed in Unit 6, 
Fourier Series. 


+ с(х)у = 0 


SAQ I 
W: page 103, Exercise 1 
(Solution on p. 23.) 

SAQ 2 

W: page 104, Exercise 10 


(Solution on p. 23.) 


10.1.2 Nonhomogeneous Problems 


The omitted passages include two complicated verifications and a reference to an 
earlier result in W which we have not asked you to read: for the purposes of this 
course the effort required to read these passages would not be worthwhile. 


Notes 


(i) W: page 127, line —9 
This means solving the differential equation by the integrating factor method. 


(i) W: page 129, Equation (29.8) 
This equation could represent the case of steady state heat conduction, with heat 
absorbed within a disc of radius R at a steady state rate F(r, 0) per unit time per 
unit area. 


(iii) W: page 129, lines 9 and 10. 
This was discussed in Section 10.1.1 (W: page 100). 


10.1.2 


(iv) W: page 130, Equations (29.9) 
о Will be asked to derive these equations in SAQ 3. It turns out however that 
Or the problem under discussion the associated homogeneous equation has 
solutions which are unbounded as гез 0. As a result the theory of Green's 


functions developed in Unit 9 must be modified by replacing the boundary 
condition 


бо = 0 
by 
G|...o is bounded. 


We shall not prove this result. 


SAQ 3 
Solve the problem 
п? 
(ran) тта = rA O<r<R, 
a,(R) = 0, 


a,(r) is bounded as r ~ 0, 
where A, is a specified function. 


HINT: Treat the cases n = 0 and n # 0 separately. The general solution to the 
associated homogeneous equation may be found in W: page 100. 


(Solution on p. 24.) 


General Comment 


In the discussion of the nonhomogeneous heat conduction problem, the first term 
on the right-hand side (W: page 127, line 6) vanishes when we put in u(0, 1) = u(r, t) = 0. 
Now suppose that the boundary conditions too are nonhomogeneous, say 


u(0, t) = f(), u(r, 0) = g(t). 


The term in square brackets would not vanish, resulting in 


n d wu c t) sin nx dx = 2[- (Nn cos ax + f()i] — n7b,(t) 
| жае sinnxdx = -[—8 ч А 


and (29.3) would be replaced by 

bio) + Put) = В) + 2 (C7 D 'ngto) nf). 
This is still a lincar first-order differential equation in one variable of the same type 
as (29.3), so that nonhomogeneous boundary conditions are still amenable to the 


finite sine transform method of solution. The solution corresponding to W: page 127, 
line —5 for nonhomogeneous boundary conditions is 


Su a Me-t) 2 "+1 ү РЕ " 
u(x.) ~ =| e^" Вт) + ш 1 * ngt) + nf (ay) dz sin nx. 
1 J0 


Example 


Solve the system 


ди д?и А 

Ae Bee х O<x<a,1>0, 

(0, 1) = 0 12 0, 

ot ls )- 120, (1) 
Ox 

u(x,0) = 0 gerrr 


Solution 


We consider first the associated homogeneous problem 


E) Pv 

S ll- O<x<ar>0, 
Q ax 

v(0, t) = t20 

дг 

aU 0) = 0 120. 

r(x, 0) = 0 0О<х<л. 


Using separation of variables, we write v(x, t) = X(x)T(t). Then 
TU). X") 


Ti^ Xe ^ ^ 99 
where 4 is a constant. The solution of 
X" - 2X = 0, 
X(0) = 0, 
Х(л) = 0, 


Х(х) = C, sin „Лх + C, cos Лх, 


with C, = O and 2 = (n + 3. (n = 0,1, 2,..) to satisfy the boundary conditions. 
Thus we obtain the eigenfunctions 
X,(x) = sin(n + 4)x forn 20. 


We can expand any piecewise smooth function g in the form* 


© 


в(х) = У d,sin(n + Dx 
nzo 


on 0 < x < n, with d, given by 


|| g(x)sin(n + 4)x dx 
d, = LL Ó 
| sin?(n + 4)x dx 
° 
agè o Г 
= Hi g(x)sin(n + 5)х dx. 
T Јо 


We can now see more clearly what type of Fourier series will enable us to tackle our 
problem (1). The associated homogeneous system suggests an expansion in terms of 
the basis (sin(n + 4)x} whatever the right-hand sides of (1) may be, provided they are 
suitably well behaved. Note, incidentally, that this is the only purpose we have in 
looking at the homogeneous system—it has no other relevance whatsoever. 


Now we can proceed as in W: Section 29 and write 
eo E 
u(x,t) ~ У b,(t)sin(n + 4)x, (2) 
120 
where 


2r 
b) = =| u(x, t)sin(n + 4)x dx, 
0 


and take the finite sine transform 
2 (*0?u : 2| du ; Ы 
ар 520% t)sin(n + Dx dx = EE: (x, t)sin(n + 3)x — u(x, O(n + Dcos(n + zi 


2p 
—(n4 4)? zf u(x, t)sin(n + 4)xdx. 
9 Jo 


* A method of proving this statement is suggested in Exercise 3 of HW’: Section 20. 
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The crucial requirement is that the term in brackets on the right-hand side can be 
evaluated, The particular choice of sin(n + 3)х, derived from the associated homo- 
Beneous system, makes this possible. 


We have, in this case 


2[(*0)u, .. Р óc 2pm 
п), axi Dsin(n + 3x dx = qn + Юл (n + pif u(x, t)sin(n + 4)x dx 
a nS „2 12 
-(-1) z^ (и + Sb, 
2 ("ди К 1Ь, 
zl gr © Osintn 4 хах = FA (0), 
and 
2[^ . 2 (-1y 
zx Meg tX | 
xL xsin(n + 5)х dx = лїї 
Thus 
db, 12 „2 1 
dU + (n + 3*5, = (—1) 41. + mpl 


which has the general solution 
= —1y" 2 I 
ba(t) = Ae Ordre { -=}14+—,5]- 
(1 +4)? n * (n + iy 


We have, from (2) together with the initial condition, 


0 = u(x,0) = y b,(O)sin(n + 4)x, 
n=0 


so that 5,(0) = 0. 


Therefore, 
(-1* 2 1 
Я = = fos 
SEES | *wegr] 


2£ (-Iy 1 ЕАР 
и) neo (n + ak tas p] MEE E 


and 


We note, of course, that the individual terms in the series are NOT solutions of the 
original partial differential equation. 


SAQ 4 
Use the finite Fourier transform method to solve the system 
2 
2 o - 5; бы) = sinx Ocxcmt»0, 
u(0, t) = u(1,1) = 0 120, 
u(x, 0) = 1 0О<х<л. 


(Solution on p. 25.) 


SAO 5 
Use the finite Fourier transform method to solve the following system (with 0 < К < 1): 
8?u Qu Pu | 
art" wat! 0«x«mt»0, 
u(0, t) = u(x, t) = 0 120, 
u(x, 0) = S4(x,0)=0 O<x<n. 


д 


This is а damped wave equation, usually called the equation of telegraphy. 


(Solution on p. 26.) 


PDE 10.1.2 


SAQ 6 


Use the finite Fourier transform method to solve 


in (0. л) x (0. я). 


y e (0. я], 
eu 
u(x, 0) = 0. ay т) = 1 x € (0. л). 


(Solution on p. 28.) 


5407 


A rudimentary, one-dimensional model of a nuclear reactor is as follows. There is a 
thin rod of active material of length a, and the reaction in it is governed by the equation 

ди _ E e 

e. 7 0c 
where u(x. t) is the neutron density and z and 3 are positive constants: the neutron 
density is taken to be zero at x = 0, x = a. The term fi?u describes the rate of genera- 
tion of fresh neutrons, and the term «?@?1/@x? describes the diffusion of neutrons. 
Show that there is a critical length of rod aa/f below which, in general, the reaction 
dies away, and above which it rapidly increases. 


+ и O<x<at>O0, 


(Solution on p. 29.) 


PDE 10.2 


Notes 


(i) W: page 132, lines —5 to —3 
You should carry this out: the new variable $ has appeared from writing 


A,(p) = Al F(p, ф)соѕ пф dp 


їг В 
вдр) = Al F(p, ó)sin nd dd 


(W: page 129, lines —9 and —8—the names of the variables must be changed, 
since r and 0 are in use already in the first equation of W: Section 30). 


(i) W: page 133, line 3 
The sense is made clearer if the words "uniformly in $" are deleted. 


(iii) W: page 133, Equation (30.4) 
Having derived this equation we use it in (30.2), putting 


"0-07 


ll 


ie. 

E 

z= 
and, also 

‚_ {r\"{0\" _ [тү 

18) \к) (А2) ° 
ie. 

gos iP. 

z= рз 


(iv) W: page 133, Equation (30.5) 
To obtain this from the previous equation, we write 
1 R 1 Ёё 1 3 
log D log(R?) 3 log(r?) 


and then combine the terms. 


The next passage in W is a verification that the expression 


R л 
u(r, 6) = Í, Г Gtr, 0:p, d)F(p, фур dp d$ 


does indeed satisfy the system (29.8). It might not have done so, since a priori there 
is no guarantee either that the system has a solution, or that such manipulations 
as the expansion of u(r, 0) as a Fourier series and the interchange of summation and 
integration leading to (30.2) are, in fact, legitimate. There is always the possibility that 
the whole operation is a purely formal exercise. The proof of validity has been omitted 
from the reading passage, since the analysis is somewhat involved, but it is important 
to realize that mechanical procedures do not always give the right answers. Here 
are two examples to underline what we mean. 
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Example 1 
Find the real roots of 
хХ+х+1=0 (1) 


(there are none!). 


“Solution™ 
Re-write the equation as 
х(х+ 0+1=0. 
But from (1) 
x*l12-x5 
so (1) becomes 
-+1=0, (2) 


which has x = 1 asa root. This, however, is not a solution of the original equation. 
[Note that x? — 1 2 (x? + x + 1)(х — 1), so solutions of (1) are also solutions of (2), 
but by increasing the degree of the polynomial we have inevitably added another 
(spurious) root.] 


Example 2 
Find all the solutions of the equation 
dy РТ: 
"m 2y xeR. 
"Solution" 


Separation of variables gives, provided y # 0 for all xe К, 


dy 
[a = fax + с, 


у= х + соогу = (х + с), 


where c is arbitrary. This represents a family of parabolas, but the left-hand half of 
any such parabola has a negative slope and satisfies 
dy 


= —2yt 
dx ib 


and so is not a solution of the original equation. Thus, we have obtained non-existent 
"solutions" as well as genuine ones. 


Thus we should realize that it is important to ask the question "Why should our 
derivation of the solution be any more reliable than these processes?" and to be 
aware that it can be answered satisfactorily. 


SAQ 8 
Show that 


222 


(5,0) > loa] re + 8 2rp cos(0 e| 


satisfies Laplace's equation for all (r, 0) and that 
(r,8) + log[r? + p? — 2rpcos(0 — ф)] 


satisfies Laplace's equation for (r, 0) # (p, ф). Hence verify the statement in W: page 
133, line —11 that (r, 0) ++ G(r, 0;p, d) satisfies Laplace's equation except at the 
point (p, $). 


(Solution on p. 30.) 
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Notes 


(i) W: page 135, lines 10 to 12 
Since D is bounded it is contained in some circle D, of radius К, say. Thus, for 
any point (r, 0) e D, 


u(r, 0) — ff, G(r, 0: p, ф)Е(р, ф}р dp аф 
D 


= [ж Lon? + p? — 2rp cos(0 — ф)] 


2,2 


+ юе + та — 2rp соѕ(0 — ero ф)р dp db 


satisfies 


LL M 1 ди _ E 


or ғ Or г? ag? 


Now in forming Vu, we can differentiate the second term in the braces under the 
integral sign; this term will then vanish, since it is harmonic on the whole of the 
region D, as we have seen in SAQ 8. Similarly, the first term vanishes under V? 
throughout that part of D not also in D. However, this term has a singularity at 
т = p, 0 = ф апа we cannot use the same procedure in D. Therefore, the function 


t 2 2 
vír, Ө) = “Tr f log[r^ + p? — 2rp cos(0 — d)|Flp, ф)р dp аф 
J D 


satisfies 


= Vu F 


throughout D. 
The condition in W that F be continuously differentiable is required so that the 
given formula for u satisfies Poisson’s equation. 
(п) W: page 136, Equation (30.11) 
This equation separates the "singular" part of the Green's function, reduced to 
its simplest form, from the smooth (harmonic) part, which is denoted by y. 
(ili) W: page 136, line — 11 
We require V? = 0. and are given the form of у on the boundary r = R. We 


know that 

r"sin n) 
and 

r"cos nO 


satisfy Laplace's equation in the circle 0 < r « R, and can therefore write down 
the form of the solution y by multiplying the nth term of the series by that multiple 
of т” which yields the correct value on the boundary. 


SAQ 9 

W: page 140, Exercise 1 
(Solution on p. 31.) 
SAQ 10 

W: page 140, Exercise 2 


(Solution on p. 32.) 


SAQ HI 
W: page 140, Exercise 5 


(Solution on p. 34.) 


10.3 GREEN'S FUNCTION: AN ALTERNATIVE APPROACH 
10.3.0 Introduction 


In this section, we seek an approach to Green's function which relates closely to 
Equations (28.4) in И: page 122. Accordingly, we return to the one-dimensional 
Green's function and show how its definition can be expressed in terms of the 
ó-function: this definition is then extended to more than one dimension. 


Finally, we apply Green's function methods to Poisson's equation in two dimensions. 


10.3.1 The Delta Function 


In Unit 9 we encountered the idea of an impulse function: somewhat analogously, in 
Unit 7, Overhead Wires we made use of a force per unit length which is zero except 
at the point of application of a force (at which point the force per unit length is 
not defined). In a similar way, we may wish to represent a point mass by a density 
function which is zero everywhere except at the mass, where it is "infinite". To permit 
a mathematical description of such physical situations, we use the symbol ó(x — č) to 
denote a quantity which is zero everywhere except at č, where it is not defined, and 
such that 


Po ap [D беба, f) 
А Гох) (х — ё) ах ZH 260 f) a) 


for any function f whose domain includes (о, f). The object ó(x — £) is known as a 
delta-function (or ó-function) and you have met it before in Unit M201 12, Linear 
Functionals: it is not really a function but a generalized function or distribution. It is 
important to realize that no meaning can be attached to ó(x — č) except in terms of 
integration. 


The idea of a delta function can be extended to higher dimensions. and we can intro- 
duce a generalized function d(r — p) in three dimensions such that 


ó(r-p)-20 ráp 


^ jf) вер 
[Jf sms - gpdV = n оёр (2) 


where / is a function whose domain includes D. (Try to construct the corresponding 
definition for the two-dimensional 5-function.) No confusion should arise out of using 
the symbol 6 for the delta function in different dimensional domains. 


It is possible to express ó(r — p)interms ofthe one-dimensional ó-function introduced 
in Equation (1). If r is expressed in Cartesian coordinates, then we replace dV in the 
integral by dx dy dz, and you can casily verify by using an arbitrary function f that 


òlr — p) = (х — &)óty = ole — eor 
where p is the vector (č, iC). Similarly, r could be expressed in spherical or cylindrical 


polar coordinates. For example, if r is expressed in spherical polar coordinates as 
(r, 0, $) and if p is (p. V, v). then we replace dV in the integral by r?sin 0 dr 40 dé. so that 


: &(r — p)d(0 — )8(ф — Dui А 
бе — р) = —— = PÜsnOÓ —— (3) 


provided that p # 0. If p = 0, we have 


б) + 


òlr) = aet 


* ]n expressions like this the equals sign means that, when placed in an arbitrary integral over a domain 
{with the correct dimension}, both sides yield the same result. 


л 


We can verify this result by looking at (2) and integrating over a sphere: 


sr Sr). 
ff m ed r^ sin dr d0 dd = ff fre. 0, on sin 0 dr d0 dd 


- mof Í ^ sin 0d0 do 
4n odo 


= f0) 


= ff f(r)ó(r) dV. 
D 
SAQ 12 


Find an expression for ó(r — p) similar to (3) in plane polar coordinates with г = (r, 0) 
and р = (p,w). What is 6(r)? 


(Solution on p. 35.) 


10.3.2 Application to Differential Equations 


To find a solution to the ordinary differential equation 

Ци)(х) = —f(x) xe (a, f) (4) 
where L is the linear differential operator given by 

L:u —»5 (ри) + qu 


we construct the Green's function G. A particular solution of (4) is now given by 
B 
19 = | fies at. 
so that 
в Hp рх (ч 
ць) = L [ fct. ad. 


Note that, since L operates on x only, it can be moved inside the integral. Comparing 
with (4), we have 


Л 
f setters 245 = -ro 


Also we know that 
ЦС](х, 2) = 0 forx Ф é, 
Thus we see that 
L[G](x, ё) = —d(x — ё). (5) 


Indeed, we can use this equation as the basis of an alternative definition of the Green’s 
function ; namely for each č e (o, fl) we specify x — G(x, č) as the continuous solution 
of Equation (5) which also satisfies the homogeneous boundary conditions 


Gla, б) = G(f, č) = 0. 
This is in fact another way of writing Equations (28.4) in W because G is a solution of 
[pG] + qG =0 for x # ё, 


and integrating (5) over the small interval (č — e, č + c) we find that 


Ete Ete 
[е | + qG = —1, 
4-г £-r 


using the integration property of the ó-function. In the limit as ¢ ~ 0. 


d 


7 ; 1 
Glo — Со = B 


in accordance with Equation (28.44) in W. Thus the 5-function provides a convenient 
notation for expressing the two conditions (28.4a) and (28.4d). 


SAQ 13 

Suppose H is given such that for each СЄ К the function 
х» Н(х, ё) xeR 

is a (continuous) solution of 
(pH'Y + qH = —ó(x — é) xeR 

and, for each xe R, 
Er H(x,& čER 

is continuous. Show that 


и) = | " Ho Os@ dé 


is a particular solution of 
(pu) + qu = —f іп Е. 
(Solution on р. 35.) 
Example 
Find G satisfying 
G'(t, x) + o?G(t,x) = —ó(t т) te R, (6) 
G(t,7) = 0 1<т, 
for each тє К. 
Hence, show that a solution of the ordinary differential equation 
и" + аи = -f ТЄР, 


is given by 


u(t) = -f JOPED a 


Solution 
The general solution of 
G'+a0°G=0 
is B sin(wt + f). We determine f from the condition G(r, т) = 0. Thus 
G(t, x) = Bsin w(t — т). 
Now, integration of the differential equation (6) from т — £to t + = yields, as c ~ 0, 


G(t + 0,7) – G(r 0,7) = —1, 


Hence 
—sin o(t — t 
—sino(t — 1) DU 


G(t,t) = а) 


0 t<t. 
Using the result of SAQ 13 a particular solution of the nonhomogencous equation is 
u(t) = f G(t, т) f(t) dt 
-a 


t sin a(t — т) 
-Í КШ б 


The generalization of these results to higher dimensions is now disarmingly simple. 
In analogy with (4). we should look at 


Ди) = —F inD. (7) 


where, to be specific, D is a two-dimensional domain, and show that 
u(r) = | G(r, p)F(p) dp reD, (8) 
D 


is a solution, where G(r, p) is a (continuous) solution of 

LIG], p) = —ó(r—- p) гєр. (9) 
(We have used dp in place of the more usual dA to stress that р, not r, is the variable 
of integration.) 
SAQ 14 
Verify that (8) is a solution of (7). 


(Solution on p. 35.) 


10.3.3 Poisson’s Equation 


We have seen in Section 10.3.2 that if 

L[G](r. p) = —d(r — p) reD (9) 
for some linear differential operator L then a solution of 

Ци) = —F in D 


is given by 
u(r) = I G(r, p)F(p) dp. (8) 
D 


Clearly, if G is chosen as that solution of (9) which also satisfies the boundary condition 
G(r, p) = 0 forrec 

where C is the boundary of D, then и (as defined above) also satisfies 
u(r) = 0 forreC. 

In this case G as defined here is just the Green's function which we met in Section 10.2. 


We show how to derive the form of G directly for the case when L is the Laplacian 
operator V^. We must solve, for each pe D, 


У2С(г, р) = —ó(r— p) гєр, (10) 
G(r.p) = 0 ` rec. 


For a given рє D we can choose a disc D, which is wholly contained in D, such that 
ped. 


We now choose the point p as the origin of a new set of polar coordinates (7, 0). 
Changing the origin does not change the Laplacian operator (this is seen best by 
looking at the definition of V? in Cartesian coordinates and noting that it is not 
altered by a change in the origin). Let 


g(t) = G(r, p) 
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where 
ro=r—p, 


ie. f gives the coordinates of r with respect to the new origin p. Then integrating over 
the disc D we have 


I, V?g(r) dA = -ff om dA = —1, 


since 0€ D. To evaluate the left-hand side of this equation we use the Divergence 
Theorem (W: page 53) for the vector field grad g: thus 


HH, Vig- n div grad g 
В б 


= g, n- grad g 
c 
where C is the circle bounding D and n is a unit outward normal. We have 
n-gradg — 0g 
cr 


since C is a circle whose centre is at the origin. 
We now require a solution of Equation (10), which must salisfy 


9g log 180g .- 


а тарар = 


subject to the additional property that 


0 for F + 0, 


f 8 _ 
LR 


Let us look first for solutions of the homogeneous differential equation which are 
independent of 0, і.е. 22/20 = 0. These may be shown to be of the form 


(7,0) = Alogr + B. 


1. (11) 


Then 
дв m onC 
OF re 


where ғо is the radius of C and 


m 2 


[ 28 Af ds= тй = 2na. 
To To Jo 


Је oF С 


From Equation (11) we require 
А = ——. 
so that 
i logr 
2л g 
is a particular solution of a 
V?g(F) = — ôf). 
Now let y be the solution of 


?у=0  inD, 


1 = 
уг: р) = x 108? rec, 


where F is the distance from p to r: then we recover the result in H^: page 136 that the 
Green's function is given by 


G(r: p) = ~ 2.180?) + (г: p). 


In polar coordinates, 

FP = r? + р? —2rpcos() — p): 
in rectangular Cartesian coordinates, 

P = (х EP + (р — 1). 
where (2, ij) gives the coordinates of p. 


If we are given a domain D which is a proper subset of R?, then we can often obtain the 
Green's function for V? on D by the following argument (much as in SAQ 9). Let P 
be the "source" point pe D. О a point not in D and T the point re D. Now, any 
multiple of log QT satisfies Laplace's equation throughout D (since Q ¢ D) and so 
does a constant. We therefore seek a suitable choice of Q and suitable values for A(p) 
and B(p) such that 


Gtr: p) = a; l-log PT + Alp)log QT + Bp] 
bI. QDE” 
7 [s ST аш 
vanishes on C. In particular we try to choose О so that QT/PT is independent of r 
when T e C. 


Example 
Find Green's function for Poisson's equation in the half plane x > 0. 


Solution 


We need G(0, y: С.р) = 0. If P is the point (2, 1), Q is (— č, 0) and Т is (x, y), then 
PT = QT when T is on the boundary x = 0, and we choose 


G(x, yin) = 
so that С = 0 on x = 0 is satisfied. Then 


1 x40 joan 
Glx, y: 6 = 7 [e E + | 


(x = 4)? + (у = п)? 


SAQ 15 (Optional) 


Let O be the centre of a circle of radius R. Let OPQ bea straight line and let Pand Q 


be such that OQ . OP = R?. Show that the ratio SP/SQ does not vary as S describes 
the circumference of the circle. 


(Solution on p. 36.) 


SAQ 16 


Use SAQ 15 to find the Green's function for Poisson's equation in the circle 
O<r<R,0<0< 2л. 


(Solution on p. 36.) 


We can use the Green's function to solve the general nonhomogeneous Poisson's 
equation as follows. We seek the solution (in two dimensions) of 


Vius-F in D, (12) 
u-f on C. 


In SAQ 11(c) of Unit 3, Elliptic and Parabolic Equations we obtained the general form 
of Green's Theorem, 


rf д ди 
23 yy? exu ees aH ge 
[Lov v — vV^u) dA $. l an v z ds, 


for any (differentiable) scalar fields u and v. Let u be the solution of our problem (12) 
and let 


v:(x, y) —9 G(X, y: Хо, Yo) 


20 
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for some fixed point (хо. yo) e D. 
Then Green's Theorem gives us 


i (uV^G — GV^u) dA = j us — e$ ds. 
D с| дп дп 


Now, we have 
V'G(p:ro) = —ó(p — гу) in D. 
where го = (хо, yo), 
Уи = —F in D, 
u=f on C 
and 
G=0 on C. 
Hence the equation above reduces to 


- ff u(p)ó(p — го) dp + ff С(р:то)Е(р) dp = £ Jt paat: 
D D JC à 


n 


и(го) = I G(to: p)F(p) dp — $ 2e (ro: Vf) ds, 
D Jc on 


using the symmetry of the Green's function (which, to be rigorous, we should have 
proved). In Cartesian coordinates we write this as 


u(x, y) = ff G(x, y: č, n)F(č, n) dé dy — 0G iy; у:ё т, nds (х, y)eD. 
D c дп 


Note that дС/дп denotes the derivative in the direction of the outward normal (on C) 
of the function 


(ё, п) > G(x, y: ё, 1). 


SAQ 17 


Obtain Poisson's integral formula (И: page 103) for the solution of Laplace's 
equation in a circle, 


Vu = 0 in the circle r < К, 
u(R, 0) = f(0) 0<0 < 21, 
using the methods of this section. 


(Solution in W: page 139, lines — 10 to —7.) 


SAQ 18 
Use the Green's function obtained in the text to show that, if 
2 2 
би бно х0, 
ôx? ду 


и(0, y) = f(y) YER, 
then a solution is given by 


ap a, 
и.) = т ах? + (у = п)? 


(Solution on p. 37.) 


PDE 10.4 
10.4 SUMMARY 


In this unit we have extended Green's functions to partial differential equations in 
two dimensions, The finite Fourier transform was introduced as a method for reducing 
a nonhomogeneous partial differential equation to a set of ordinary differential equa- 
tions, and shown to lead to Green's function for Poisson's equation. 


We then demonstrated how Green's function could be defined using ó-functions. 
Poisson's equation in two dimensions was treated in this way with nonhomogeneous 
boundary conditions, making use of the Divergence Theorem and Green's Theorem. 
The ó-function was defined by 


òr- p) 20 г р, 
Е | Уе рер, 
„1980 = p)dr = | péD, 


for all functions f. (This definition may be interpreted in one or more dimensions.) 
We saw that the Green's function for the problem 
Lu = -F in D 
is given by the continuous solution of 
LG(r:p) = —d(r — p) reD 
G(r: p) = 0 rec 


for each pe D, where C is the boundary of D. 
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10.5 SOLUTIONS TO SELF-ASSESSMENT QUESTIONS 


Solution to SAQ 1 


We need to express sin? in the form 


» 
зао + У, (a,cos nO + b,sin n0). 
n=l 


Now, 
isin 0 = de — e7"), 
so, taking the cube of each side, we obtain 


—isin?0 = Her” — Зе? 4 3970 — e730), 


sinü = — l2 sin 30 — 6 sin 0) 
= {зіп 0 — ksin 30. l (1) 


For ғ < 1, the solution is given by Equation (24.6) in W: page 101, 
u(r, 0) = 3ag + È rla,cos nÜ + b,sin n0), 

with the coefficients given by 
u(1,0) = зао + У tancos nO + b,sin n0). 


Comparing with (1), we see that 


u(r, 0) = 13r sin 0 — r?sin 30) йк 1. 


Solution to SAQ 2 


To solve 
2 2 
уе ltp pI ug 1<г< 8,0 < 0 < 2л. 
д rór 1? 00° 
(1,0) = 0, 
u(R, 0) = f(0). 
we рш u(r, 0) = w(r)O(0). 
Then the equation becomes 
"dla ЧГ P 
rw" rw(o)  — O"(0) _ п?, say. 


Ww т w(r) ^ e(0) 
We now have raw” + rw’ — n?w = 0, whose general solution is 
wol) = ao + bolog r п = 0, 
walt) = a,r" + br" n# 0, 
and Ө” + п20 = 0, with the eigenfunctions (for solutions with period 2л) 
1.0 n = 0), ү 
cos пф, sin nO neZ*. 
The boundary condition u(1, 0) = 0 is equivalent to w(1) = 0, so that 
wo(r) = bolog r. 
wur) = br" т) neZ*. 
Hence we seek a series solution in the form 


u(r.0) = bologr + Y. (r^" — r')(c,cos nO + d,sin n0). 
n= 1 


n= 
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To satisfy the boundary condition at r 2 R we require that 
Д0) = bolog R + У (К " — R")(c,cos nO + d,sin nO) 0 € [0, 2л). 


wel 
Thus 
bolog R = 1A; 
(R^" — К"), = A, nez* 
(R^" — R"\d, = B, neZ*, 


where А,, B, are the Fourier coefficients 


A, = if Гф)соѕ пф аф п 2 0, 
T-n 


В, = HN nó аф п> 0. 


We consider the function 


© 


logr D 
u(r, 0) = 3А ER + 2 RIR 


(A,cos n0 + B,sin n0) 


1<г<К,0<0<2т. 
If 


1f 
c=i f voa, 
7-5 
then |A4,| < С, |В, < С, and the series for и and its first and second partial derivatives 
are dominated by 


for 1 <r € rg where rg < R. Thus, by Weierstrass’ M-test the series for u(r, 0) 
and its first and second partial derivatives are uniformly convergent for 1 <r < го. 
Hence we may differentiate the series term by term and verify that и satisfies Laplace's 
equation for 1 < г < К, and also и is continuous for 1 < r < R. By an argument analo- 
gous to that given in W: page 101, lines — 13 to —4 we may show that it is in fact 
continuous for 1 < r < К and so the sum of the series given above for u(r, 0) is the 
solution to our problem. 


Solution to SAQ 3 


To solve the problem we construct the Green's function G, which must satisfy the 
following equations for each рє (0, R): 


n? 
(Gy —- —6=0 r$ р, 1) 
С, o is bounded, С|,., = 0, 2) 
G(r, p) = Ср, р), (3) 
7 i 1 
G оо Е б шев = Up (4) 


Note that we have replaced the continuity condition by the symmetry condition as 
the latter is easier to apply. 


For the casen = 0а linearly independent pair of solutions of(1) is (1, log r). To satisfy 
(2) we set 
С,(0) (logr — log R) r 


би») = { и 
ОР с) r<p 


24 
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Condition (3) now gives us 
Cilp) = C, 
C3(p) = C(log p — log R), 

where C is some constant which we determine from (4), obtaining 
Çai 

Thus the Green’s function is 


log R — logr r 


2 
боб", p) = { 2p 
olr, p) log R — log p rsp. 


and the solution to the problem is 


R 
aor) = |. баў ple Aotpdp 


"К R R 
= | log 7 pAo(p) dp + Í log 7 pAd(p) dp Oxr«R. 


For пФ Owe may assume without loss of generality that n > 0.A linearly independent 
pair of solutions is r^, т^", and to satisfy the boundary conditions (2) we require 


ale ste 
ato Ie 0] >» 


Рр)" r & p. 
Condition (3) now gives us 
Di(p) = рр", 


СЕТЕ 


and we determine the constant D from (4), obtaining 


M 
~ 2nR" 


Thus the Green’s function is 

оү {кү fr | 
2n\R} | ir. R 

1 r n R n p m" Е 
so - e] "s^ 


and the solution to the problem is 


Gr, р) = 


R 
a,(r) = i Сы, DPA(9) dp 
0 


cai LE - Ge) ЕЛ once] 


О< г< К. 


Solution to SAQ 4 " 


Separation of variables for the corresponding homogeneous system leads to the 
boundary value problem 


X"--AX = 0 in(0,z) 
X(0) = X(x) = 0, 
which provides a suitable basis for expansion: the eigenfunctions are 
{sin nx, n = 1,2,...}. 


Therefore we write 
© 

ибх, 0) ~ У b,(t)sin nx, 
1 
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where 
л 
bt) = 3i u(x, t)sin nx dx. 
T Jo 


The finite Fourier transforms are given by 


"Ou. 2[ ôu . GNE FX ee 
= | zasinmxdx —-|z-sinnx — uncosnx | ——n u sin nx dx 
g 0x л| Ox о п o 
= —n’b,(t), using the boundary conditions; 
2 "ди. 4, . abt), 
al at sin nx dx = dO 
and 


2 [^ nrd dx = | qd 
л Jo : FUR 0 mex 


(sin x is its own Fourier Series, consisting of one term only.) Taking the finite sine 
transform of the differential equation we have 


db, „ [i n=, 
di * b 00 nz2. 


dt 
The initial condition, with t = 0, gives 
u(x, 0) = 1 ~ У b,(0)sin nx. 
l 


(1) 


The Fourier coefficients are given by 


0 neven, 
2". 
ыб) = = | 1. sin me dx = Lm 
пл 


and, in conjunction with (1), we obtain 
T e -t 
bi) = =е '-(1—e79, 
bt) = 0 neven, 
LA 
b(t) = —e^"* nodd,n» 1. 
пт 
Тһеп 


єй 4 ү 
u(x, t) = (1 — e~sinx + - Y, —е "sin nx. 
noda P 


Solution to SAQ 5 


As in the solution to SAQ 4, we expect to be able to use the expansion 


u(x,t) ~ Y. b,(t)sin nx. 


п=1 


th 


xy 
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Accordingly, using the boundary c itions at x = у! ai n 
a onditions at x finite si 
i y 0, п, we obtain the finite sine 


2 [78u d 2 
=] scysinnxdx = — 
п), Ox? x n*b,, 


27 Ou. d?b 
= f gz sin пх dx = dg 


2[*Qu. db, 
al gr Sin nx dx =з" 


0 neven, 


2T. 

zf зіп nxdx = 4 4 

Tjo — nodd. 
пт 


Therefore, multiplying the differential equation by = sin nx and integrating, we obtain 
= А 


itp " 0 neven, 
"+ 2k- + nb, = 4 
а? а" m” odd. 


The general solution is easily found to be 


0 n even, 
bt) = e™ [A cos /(n? — K?)t + Bsin /(i* — Kr] + 
XC odd. 


The initial conditions yield 
db, 
5,0) = 70) = 0 


whence 


0 neven, 
0= А4, + 4 
p n odd, 


0 = —kA, + /(? — К?)В,. 
Thus 

А, = В, = 0 п even, 
4 —4k 
3 


mz Pe 3 taht 
n n Jn k?) 


n odd, 


and, finally 


u(x, )= 4 У пеле = к); ра" - ex| - i sin nx. 


nodd 
This problem represents a form of damped vibrations of a string driven by a uni- 
formly distributed force along its length and starting from rest. Physically we would 
expect the vibrations to damp out and the string to settle into an equilibrium position 
as t ~~ 20. Thus ? 


" 4 JM 

lim u(x, t) = = У, -sin nx 

ге T nodd 
describes the “final” shape of the string. Since this is independent of t, we should 
also expect it to result from writing ĝu/ðt = 0 in the original equation. Let 


их e lim u(x,t): 


tna 


then u, should satisfy 


2 
d'uy 
dx? 


u,(0) = и„(л) = 0. 


Hence 


and the boundary conditions give A = іл, В = 0, so that 
u(x) = —ix(x — n). 
We expand this as a sine series over [0, л]: 


T 
и„(х) ~ У D,sin nx, 
п= 1 


where 
2 [* К 2 7 
D, == | —4x(x — n)sin nx dx = — | -= cos nx 
ajo ^ nn 0 
0 n even, 
=) 4 
— nodd. 
mn 
Hence, 
U(X) = —ix(x — л) = Y, sin nx, 


as expected. 


Solution to SAQ 6 


PDE 10 SAQ 5-6 


As in SAQs 4 and 5 we take (sin nx} as a basis. If и is a solution we may expand it 


in a Fourier series as 


u(x, y) ~ У b,(y)sin пх, 


where 
2[" : 
Ьу) = zf u(x, y)sin nx dx. 
Tjo 
The differential equation becomes, on transformation, 
dh, db, 
dy n^b, + b, = dy! 


the boundary conditions become 


- (? — 1), = 0; 


u(x,0) = 0 = Y. b,(O)sin nx, 


п=1 
ai со 
Pree n=ls 2. bi(n)sin nx, 


so that 
5,(0) = 0 all n, 


0 n even, 
п) -| 4 


пл n odd. 
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y^ 


The general solution of (1) is 

by) = AP Вее 
forn = 2,3,4..., whilst 

by(y) = А, + Byy- 
The boundary conditions give us 


A, + В, = 0 ж==& E .„, 
0 п еуеп 


(n? — 1) [A e — Be] д 
sn n-—3,5,7,...5 


in addition, 
and 


Thus 


A, = В, = 0 neven 


2 


~ an(n? — 1) cosh x(n? — 1) „ш 


The solution is therefore given by 


rs 


a. 4 1 sinh y(n? — D* , 
u(x, y) = 2) sin x + RN ie — TF cosh nin? — D sin nx, 


пода 


since the series may be shown to be uniformly convergent on [0, л] x [0, я]. 


Solution to SAQ 7 
The basis for our expansion will be the eigenfunctions of the system 
РХ 
2 + AX = 0 x e (0, a), 
dx 


X(0) = X(a) = 0, 


obtained by separating variables in the given equation. The functions 
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are suitable, where, to satisfy X(a) — 0, we require 


that is, 


Hence we put 
= . ATX 
u(x,t) ~ Y, b(t)sin — 
п=1 а 
with 
2[* . ATX 
bt) =- [ u(x, t)sin —— dx. 
a Jo a 
Then, 
2[^0?u . nnx nên? 
= | —>sin— dx = ————Ь„, 
a Jo Ox? a а?" 
2 "ĝu . тх db 
=| —sin—dx =—- 
а), д а dt 
The equation becomes, under transformation, 


db п2д?о? 
~=- b, 2p. 
qi Th, Bb, 


or 


for each ne Z*, which has the general solution 


1 


The coefficients B, are determined by the initial conditions. In general, B, + 0 and 
it is immediately apparent that, if 


I?n?o? 
- ; > 0, 


nna? 


b) = вело [i = 


а? 


p 


a 
the solution tends to infinity. (Note that as t increases, the first term dominates the 
series.) Thus the condition for the reaction to grow is 
S 
a Г ` 
since a, B > 0. 


In physical terms, the coefficient f, which governs the rate of production of fresh 
neutrons, predominates over «,.which governs the rate of diffusion out of the bar at 
its ends. 


Solution to SAQ 8 
We shall show that 
h:(r,8) — Іор [ f(x, 0)] 
is a solution of Laplace's equation wherever f(r, 0) + 0, for the functions 
rp? 


fi: 0) — R? + тт ce cos(0 — ф) 


30 


and 


Л 0) — r? + p? — 2rpcos(0 — ф), 


where p and ф are constants. 


We have 

ah lof 

or fôr 

2h 127 1[8fV 

а? for EA 
and 

Ph 1 Of 1 әд 

200 f OOF л 
Непсе 

2 

Е 

Now 


V?[rcos(0 — ф)] = 0, 
as we saw in Section 10.1.1. Let 

f:(r, 0) э P + Qr? — "rpcos(0 — ф), 
where P, Q are constants. Then 


4[Qr — p cos(0 — ф)]? + 4p?sin*(0 — ф) 


V?h(r,0) 


4Q 
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f(r, 8) 
_ A[PQ + Q^r? — 20rp cos(0 — Ф) — 02° + 20гр cos(8 — $) — p°] 


Lir, O° 


_ PQ — р?], 


Ut, 


oy 


Ut, 9)? 


For f, we have Р = R? and Q = p?/R?; for f; we have P = p? and Q = 1. In either 
case PQ — p? and so 


Vh = 0. 


Thus 


V2G(r, Ө; p, $) = Vlog fir, 0)  V?log/a(r, 8) 


Solution to SAQ 9 
We require a Green’s function G such that, for each (p, d) €(0,R) x (0,7), 


-0 


for (r, 6) + (p, Ф). 


G(r,0: p, $) = -Š logt? + p? — 2rpcos(0 — ф)] + y(r, 8: p, Ф), 


where 


У2у = 0 


and 


(r, 0) e (0, R) x (0, x) 


G(R, 0: p, $) = G(r,0:p, à) = Gr, n: p, d) = 0. 
Let G denote Green's function for the circle 0 < r < R given by Equation (30.5) in W: 


1 
G(r, 0; p, Ф) = а= 


4n 


[mte + p? — 21р соѕ(0 — ф)] 


ee + 


729? 


7 


2,2 
RO 2rp cos(0 — | | 
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Then 
V?G(r.0: p.) = 0 (0.0) (p. d) 
and 
V?G(0:p.—d)20 (7.0) (p, — ф). 


Thus the sum of G(r, 0: p. ф) and any multiple of G(r, 0: p, — ф) satisfies Laplace's 
equation at all points of the semicircle 0 < r < К, 0 < 0 < x except (r, 0) = (p, $). 
and vanishes for ғ = А. We seek a suitable 2 such that 


Gir, 0: p, $) + 2G(r,0: p, 0) = 0 0 = 0, я. 


Thus we require 


cos(—$) = —2cos $ 
and 
соѕ(л — ф) = —2cos(n + ф), 
which are satisfied by 4 = — 1. Hence Green's function for the semicircle is given by 


2 


uc ad 
G(r, 0: p, ф) = x t log(r? + p? — 2rpcos(0 — ф)]+ log [e + a — 2rp соѕ(0 — ^| 


+ log[r? + p? — 2rp соз(0+ ф)] — led ie +E — arp cos(0 + a} 4 


which is in the required form, by SAQ 8. 


Solution to SAQ 10 


We may proceed as іп W: Section 29, to obtain the equations at the top of W: page 130 
with the additional boundary conditions 


a,(1) = b,(1) = 0. 


We must then construct the Green’s function g, given by 


ү п? 
(ren —-=0 гр, 
8.1 = КАК = 0, 
grs p) = gp. г), 

1 


p 


These equations may be solved in a manner analogous to that used in SAQ 3. 


Ele oso — Bn b=p-0 = 


When п = 0, the general solution of the homogeneous equation is 
Alogr + B, 


50 


Ai(p)log ғ <р, 
Boll’, p) = R й 
Ayplogz г> р. 


The symmetry condition gives 
R 
A,(p) = A log p 


Alp) = A log p, 
and, since 
I 


go|r-5*o — £ol-25-o0 = — — 
p 


a 
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we have 
Alogp A log(R/p) _ 1 Я 
p р р n log R^ 
so that 
log r log(R/p) 
log R 


Bol", p) = m 


log R 


For n > 0, we obtain, in the same way 
Dp RY py . 
IR" R| (5) TIR r Xp, 
р" рт" RW py | 
ишек [| -(}] т>» 


Arguing as in W: Section 30 we find that 


810% р) = 


1 
G(r, 0: p, Ф) = = Bolt, p) + у Letro )cos п(0 — ф) 


1 flogrlosiR/p), 2 тс" TIRY [py | 
x log А P. E EEE [I аа-а) Pes 


lflogplogR/], 2 р" p~” RV ew | 
A Tog R tà AO x) | х (s) Je n(0 — ф) ) rp. 


An alternative method of solution is to proceed as in W: page 136. We write 


С=у*+у 


where 


1 2442 
7%(,0;р,ф) = — a log[p? + r? — 2pr cos(0 — 4)] 


1 э py)" 
= el t - У 4 cos n(0 — | r<p, 


n= 


z i T» Y: е)" cos п(0 — | rp. 
n= 1! 
To find 7 we solve, by separation of variables, the problem 
V?y=0 —intheannulusl <r < R, 
1 v 1 74 PES А 
70,0:р, ф) = "е - У "L cos n(Ü 9| 


a=1 


ARO: = | OB b 318) cos n(6 a): 


= 1 
760,0: p, b) = (r,0 + 2n: р, ф). 


We then obtain the series for G given above. 
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Solution to SAQ 11 


Proceeding as in SAQ 10, we look at the system 


ne п? 
(„у —-у&=0 rA 


£l, ~o is bounded, g,’ + g|.-a = 0, 


Enl", р) = gp. r), 
7 " 1 
En b=e+0 — £n [= = E 
When n = 0, the general solution of the differential equation is 


Alogr + B, 
so to satisfy the boundary conditions we require 


Ар) r<p, 


go( p) = мов) zi Box 


The symmetry condition gives 


= p| 1 
asalet] 
Ар) = А, 
and we determine the constant A by the discontinuity in the derivative, 


A 1 wok, ie A= —1. 
р р 


For п > 0, the general solution is 
Ar + Вт" 
so to satisfy the boundary conditions we require 


Alp)" r «p, 


80% р) = (o) R|[r)" | [R|" m 
Bip exl + “2 


Also the symmetry condition.gives 


2 n—R\{p\" RY MU 
A,(p) = di TE (5 + В | B,(p) = Ср", 


and the jump condition gives 


{ке е An] - Lr] + Be} -- 


1 п R 
ama" ~ В) (a + (п + ШЕ | 19 кєр; 


экю e-o e mE] , 


Following the argument of W: Section 30, we have 


so that 


gr, p) = 
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1,0: p,@)= a Bol p)+ Li ET , p)cos n(0 — d) 


1 


Solution to SAQ 12 


The element of integration in plane polar coordinates is dA = r dr d0, so that 


spa io = =: 
For p = 0, 


900). 


2nr 


б(т) = 


For, integrating over a domain D which contains the origin, we have 


If, ftr, 020, a= rof 40 


= f(0) = I ftt) dA 


which is the correct result. 


Solution to SAQ 13 
Let 


w= | 7 Hos Ofte) dt, 


and write L for the linear operator 
L:u > (риу + qu, 
so that 


Lb) = L | " HOSO ade 


| цн ӘӘ 


-f х — ƏS dé 


= —f\(x) xeR. 


{1 ROT. 1 r\" 
E TIRT dn Pr (n юа 2 ++? ЩЕ cosn — 9) 


1 R e 1 n "j n 
x + log z +5, ee AR С - viz +(n+R) |: | cos n(@ — of 


Thus, irrespective of the boundary conditions on H, we obtain a particular solution 


for u. (Note that this justifies our argument in note (iv) of Section 10.1.2.) 


Solution to SAQ 14 


u(t) = Í б\г: pfp) dp 
D 


so that 
1(и)р = LÍ G(r: p)f(p) dp 
D 


| LIG): pf (9) dp 
D 


Ш 


= f б\т — р)/ (о) dp 
JD 


=) гер. 
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г< р, 


г> р. 


Solution to SAQ 15 
Using the cosine law. where OP = p. OQ = q. we have 
(SP)? = R? + p? — 2Rp cos В. 
(SQ)? = R? + g? — 2Rq cos f. 
$ 


БА 


Eliminating cos fi we obtain 
2 R? SP 2 R? 
(SQ | TE Ba ZI 


q q p р 
Now pq — R?, so substituting for R?, 
(SQ)? (SP? 
"E (p+4q)= п (9 + р). 
Непсе, 
SP _ e =? 
SQ q R 


which is independent of ff and SP/SQ does not vary as S describes the circumference 
of the circle. 


Solution to SAQ 16 


Take P to be the point (p, ф), so that О is (R?/p, ф), and take T to be (r, 0), and, by 
SAQ 15, 


QT -R 
PT p 
when T is on the boundary. Then 


1 T R 
Gir, 0: p, ф) = 2 - a | 


satisfies the equation 


8G 18G 1 &G 
y + т ot тї 90 = ~ or 0) 


and vanishes on r = R. Therefore 


Gr, 0:p, d) = a {toe As (om? ы ату} 


me 04 


p 
R? 


р 


1 2 1 
== log(r? + p? — 2rp соѕ(0 — ф)] + — os| R + 


= — 2rp cos(0 — e| 


as found in W: Section 30. 


Solution to SAQ 17 
See W: page 139, lines — 10 to —7. 
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Solution to SAQ 18 


PDE 10 SAQ 18 


The boundary of the domain is x — 0, so we have, using the formula given in the text, 


дп 


* 0G m 
u(x, у) = -Í => (x, у: e pfo dn. 


Remember we need the directional derivative along the outward normal, so 


әс 8G 


ón дё 
= nmi Reb ык 
пр + (y = 2] 
Then, 


ЕЕ T Jn du 


-elx* + (y 007] 


__ 4 lio 
ear дё Ал 


(x + & + (у = п)? 
Sle Fy | eo 
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conducting hodies, and is an immediate consequence of what has preceded. For leta, y, z, 
be the reciangular co-ordinates of any particle р in the interior of one of the bodies; then. 


wi (2) 


be the force with which p is impelled in the direction of the co-ordinate x, 


1 ; ae D n 
and tending to increase it. In the same way — 92 and — £% will be the forces in y aud 


2 
a, and since the fluid is in equilibrium all these forces are equal to sero» hence 
Wap Vy Wa 
Gaiety 4, 


which equation Being integrated gives 
сове. 

‘This value of V being substituted in the equation (1) of the preceding number gives 
: o=0, 
and consequently shows, that the density of the electricity at any point in the interior of 
any body in the system is equal to aero. ` 

The same equation (1) will give the value of » the density of the electricity in the 
interior of any of the bodies, when there are not perfect conductors, provided we can 
ascertain the value of the potential function J^ in their interior. 


(8) Before proceeding to make known some relations which exist between the density 
of the electric fluid «t the surfaces of bodies, and the corresponding values of the potential 
functions within and without those surfaces, the electric fluid being confined to them 
alone, we shall in the first place, lay down a general theorem which will afterwards be 
very useful tous. This theorem may be thus enunciated: 

_ Let U and V he two continuous functions of the rectangular co-ordinates a, y, z, whose 
differential co-efficients do not become infinite at any point within a solid body of any 
form whatever; then will | 


Semra Го 2) ftir fi der (GE ; 


the triple integrals extending over the whole interior of thc body, and those relative to de, 

over its surface, of which d» represents an element: de being an infinitely small line per. 

pendicular to the surface, and measured from this surface towards the interior of the body. 
To prove this let us consider the triple integral à 


Je OC CCS GE}. 


‘The method of integration by parts, reduces this to 
ал du л 
J'en Sier E а аи T ferg 


dU p., ay FU dU ас 
fiam t2 fe. d Semir E EL b 
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эзох, in the commencement of his first memoir (Mem. de Institut. 1811), has inciden- 
AE given а method for determining the distribution of electricity on the surface of a 
spheroid of any form, which would naturally present itself to a person occupied in these 
researches, being in fact nothing more than the ordinary one noticed in our iutroductory 
observations, as requiring the resolution of the equation (a). Instead however of "p 
posing, as we have done, that the point p must be upon the surface, in order that the 
equation may subsist, M. Porssow dede himself of a general fact, which was then sup- 
ported by experiment only, has conceived the equation to hold good wherever this point 
may be situated, provided it is within the spheroid, but even with this extension the me- 
thod is liable to the same objection as before- - E 

Considering how desirable it was that a power of universal agency, like electricity, 
should, as far аз possible, be submitted to calculation, and reflecting on the advantages that 
arise in the solution of many difficult problems, from dispensing altogether with a par- 
ticular examination of each of the forces which actuate the various bodies in any system, 
by confining the attention solely to that pecoliar function on whose differentials they all 
depend, I was induced to try whether it would be possible to discover any general relations, 
existing between this function and the quantities of electricity in the bodies producing it. 
The М 'LAPLACE had derived in the third book of the Mecanique Celeste, from 
the use of a partial differential equation of the secoud order, there given, were too marked 
to escape the notice of any one engaged with the present subject, and naturally served to 
suggest that this equation might be made subservient to the object I had in vi ecol- 
lecting, after some attempts to accomplish it, that previous researches on partial differential 
equations, had shown me the necessity of attending to what have, in this Essay, been 
denominated the singular values of functions, I found, by combining this consideration 
with the preceding, that the resulting method was capable of being applied with great 
advantage to the Fecia theory, and was thus, in a short time, enabled to demonstrate 
the general formule contained ia the preliminary part of the Essay. The remaining part 
caret М E. principally as farnishing particular examples of the use of these 
general formulæ; their number might with great ease have been increased, but those 
which are given, it is hoped, will snffice to point ont to mathematicians, the mode of apply- 
ing the preliminary results to any case they may wish to investigate. The hypotheses on 
which the received theory of magnetism is founded, are by no means so certain as the facts 
on which the electrical theory reste; it is however not the less necessary to have the means 
of submitting them to calculation, for the only way that appears open to us in the investi- 
gation of these subjects, which seem as it were desirous to conceal themselves from our 
View, is to form the most probable hypotheses we can, to deduce rigorously the conse- 
quences which flow from them, and to examine whether such consequences agree numeri- 
cally with accurate experiments. 

е applications of analysis to the physical Sciences, have the double advantage of 
manifesting the extraordinary powers of this wonderfal instrument of thought, aud at the 
same time of serving to increase them; numberless are the instances of the truth of this 
assertion. То select one we may remark, that M. Fourrer, by his investigations relative 
to heat, has not only discovered the general eqnations on which its motion depends, but 
has likewise been led to new analytical formule, by whose aid M. M. Caucay & Potsson 
have been enabled to give the complete theory of the motion of the waves in an indefinitely 
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the accents over the quantities indicating, as usual, the values of those quantities at the 
limits of the integral, which in the present case are on the surface of the body, over whose 
interior the triple integrals are supposed to extend. 


д 
Let us now consider the part Гб “Z duo to the greater values of а. It is casy to 


see since dw is every where perpendicular to the surface of the solid, that if de? be the 
clement of this surface corresponding to dude, we shall have 


pu 


du _ dr ar 
Jm ea fuent. 


" prs p gy 
In like manner it js seen, that in the part fa Фраер" ^ due to the smaller values of z, 


and hence by substitution 


we shall have dyie= +" dy, und consequently 


au de at! 
finer ae Je a E 
"Then, since the sum of the elements represented by dz; together with those represented 
by ae", constitute the whole surface of the body, we have by adding these two parts 


тага 


where the integral relative to de is supposed to extend over the whole surface, and de to 
be the increment of т corresponding to the increment de. 


In precisely the same way we have 
к npe fe 


and 1 fisa ( v га = 


therefore, the sum of all the double integrals in the expression before given will be ob 
tained by adding together the three parts just found; we shall thus have с 


dU de | dU d, | dU d dU 
Je de de t dy de * di a= fore, 
aU 
where V and да represent the values at the surface of the body. Hence, the integral 
dV dU dV dU dV av) 
Jine a ae aE d 

by using the characteristic 3 in order to abridge the expression, becomes 
РЕ РА 


vni 

extended fluid. ‘he same formule havo also put us in possession of the solutions of 
many other interesting problems too numerous to be detailed ћего. It must certainly be 
regarded as a pleasing prospect to analists, that at а time when astronomy, from the state 
of perfection to which it has attained, leaves little room for farther applications of their 
art, the rest of the physical sciences should show themselves daily more and more willing 
to submit to it; and, amongst other things, probably the theory that supposes light to de- 
pend on the undulations of a luminiferous Suid, and to which the celebrated Dr. T. Youre 
has given such plausibility, may furnish a useful subject of research, by affording new 
opportunities of applying the general theory of the motion of fluids. The number of 
these opportunities can scarcely be too great, as it must be evident to those who have 
examined the subject, that, although we have long been in possession of the general equa- 
tions on which this kind of motion depends, we are not yet well acquainted with the 
various limitations it will be necessary to introduce, in order to adapt them to the different 
physical circumstances which may occur. 

Shonld the present Essay tend in any way to facilitate the application of analysis to опе 
of the most interesting of the physical sciences, the author will deem himself amply repaid 
for any labour he may have bestowed проп it; and it is hoped the difficulty of the subject 
will incline mathematicians to read this work with indulgence, more particularly when 
they are informed that it was written by а young man, who has been obliged to obtain the 
little knowledge he possesses, at such intérvals and by such means, as other indispensable. 
avocations which offer but few opportunities of mental improvement, afforded. 
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